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CHAPTER  I 


INTRODUCTION 

Several  Radar  Target  IdentiScaiion  (RTI)  techniques  have  b^n  developed  at 
The  Ohio  State  Uuivernty  in  recent  years.  Using  the  HectroS-ien<e  Laboratory 
compact  range  a  large  data  bas:  of  a>hCTait  RCS  measumi^ts  has  been  con¬ 
structed  for  ^verai  typo  of  targets  (airerail,  ships,  uid  ground  vehicles)  at  a 
variety  of  polarizations,  aspect  aa^o,  and  frequency  bands.  This  extendve  data 
base  has  been  used  to  analy^  the  performance  of  ^veral  differmit  clasdScation 
algorithms  through  the  use  of  computer  Emulations. 

In  order  to  optimize  classification  performance,  Ksienski  |1]  concluded  the 
radar  fiequency  range  should  lie  in  the  Hayidgh-resoaance  frequency  range,  where 
the  wavdrogth  is  on  the  order  of  or  larger  than  the  target  sze.  For  aircraft 
and  ships  with  general  dimension?  sc  the  cjrder  of  10  meters  to  KM)  meters  it 
is  apparent  that  the  High  Frequency  (HP)  band  provides  optimal  classification 
perfonaaace.  Since  existing  HF  radars  ere  currently  being  n%d  for  detection  and 
tracking  of  mreraft  and  ships  of  th^^  dimendems,  it  is  natural  to  further  investigate 
the  posdbility  of  udng  these  Gdsliug  rt^ars  as  the  measurement  devices  in  a  radar 
target  clasdfication  system. 

-4  general  radar  target  dasdficaiion  system  is  shown  in  Figure  1.  The  external 
environment  coiidsts  of  a  taiget  widcli  is  to  be  ciasdfied,  clutter  source{s),  noise 
source(s),  and  possible  calibration  refer^nce/s).  The  environment  is  observed  by 


the  multifrequency  radar  through  the  propagation  media.  The  signal  processor 
provides  radar  data  which,  when  combined  with  a  priori  knowledge  of  the  external 
environment  and  knowledge  of  the  radar  system  and  propagation  media,  allows  the 
estimation  of  target  features  to  be  used  for  classification.  Th^e  estimated  features 
of  the  unknown  target  are  compared  with  high  accuracy  features  of  sever:**  *5-igcts 
in  a  catalog  set.  The  classification  algorithm  provides  a  measure  of  the  similarity 
of  the  unknown  target  to  each  catalog  target.  The  catalc^  target  whose  features 
most  closely  repr^ent  that  of  the  unknown  targets’  feature  is  chosen  as  ihe  type 
of  the  unknown  target. 

This  study  investigated  various  propagation  conditions  which  allow  different 
target  features  to  be  estimated.  These  conditions  were  syntheszed  by  various 
cii2innel  models.  The  incorporation  of  these  channd  models  into  a  radar  target 
identification  system  computer  simulation  provided  estimate  of  resulting  deifi¬ 
cation  performance.  As  a  result,  the  relationship  between  channel  conditions  and 
classification  performance  were  found. 

HF  radars  are  generally  categorized  by  the  propa^tlon  mods  by  wiueb  th^ 
are  designed  to  operate,  either  1)  surface  wave  or  2)  sk3rwave.  Each  propagation 
mode  incites  different  characteristics  into  the  dassification  system,  and  th»efore 
should  be  investigated  independently. 

Because  of  their  unique  propagation  medianism,  surface  wave  radms  have  the 
ability  to  look  beyond  the  horizon  with  detection  ranges  op  to  300  km.  Surface 
wave  attenuation  is  a  function  of  many  paraiiieters,  most  impottantly  ground  con¬ 
duction  and  freq-'ency  [6].  !n  order  to  achieve  practical  detection  ranges  for  given 
transmitted  power  levels  the  radar  wave  usually  must  propagate  over  water.  Fur¬ 
thermore  only  a  vertically  polarized  wave  propagates  significantij,  thus  limiting 
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Figure  1:  Gci,cral  classification  svsteni. 


the  polarization  diversity  of  our  measurement  system. 

HF  skywave  radar  utilizes  ionospheric  reflection  as  the  propagation  mecha¬ 
nism.  The  skywave  radar  has  an  average  detection  range  of  1000  km  to  4000  km, 
which  when  coupled  with  an  antenna  steerable  in  azimuth  up  to  120°,  provides 
coverage  of  over  15  million  square  kilometers  [12].  The  ionosphere  supports  HF 
radio  propagation  at  altitudes  of  100  to  500  km  and  consists  of  i.>nized  regions 
referred  to  as  kyers  which  are  commonly  labeled  D,  E,  Fi,  and  F2.  Ions  at  these 
altitudes  are  mainly  produced  by  solar  radiation  which  includes  particle  radiation, 
ultraviolet  light,  and  x-rays.  In  general,  it  is  the  electron  density  distribution  at  a 
given  altitude  which  effects  the  propagation  of  the  EM  wave  to  the  greatest  extent. 

The  actual  ion  distribution  at  any  time  is  a  function  of  many  factors  includ¬ 
ing  source  radiation  intensity,  angle  between  zenith  and  solar  radiation  rays,  ion 
distributions,  and  various  ion  combination  mechanisms  as  described  in  Davies  [14]. 
Therefore  the  ionosphere  is  continually  changing  and  must  be  monitored  in  real¬ 
time  to  accurately  determine  its  present  structure.  It  is  the  variation  in  electron 
density  with  altitude  which  causes  variation  in  the  refractive  index  and  results 
in  the  bending  of  the  ray  path  and  its  return  to  earth.  It  is  important  to  note 
the  ionosphere  only  supports  propagation  over  a  limited  band  of  frequencies,  and 
signals  of  different  frequencies  may  travel  along  different  ray  paths. 

The  earth’s  magnetic  field  also  affects  the  propagation  of  the  radar  wave.  The 
magnetic  field  and  ionosphere  produce  a  magnetoionic  medium  which  supports 
propagation  modes  with  spedfic  polarizations.  These  polarizations  are  a  function 
of  the  magnetic  field,  parameters  of  the  ionosphere,  and  the  direction  of  propaga¬ 
tion  as  given  in  Davies  [14].  As  a  result,  a  transmitted  wave  of  a  given  polarization 
may  split  into  two  ray  paths  which  may  or  may  not  travel  the  same  route  to 
the  target.  Therefore,  in  general,  the  polarization  of  the  transmitted  wave  is  not 


necessarily  the  polarization  of  the  wave  v/hich  illuminates  the  target,  and  the  po¬ 
larization  of  the  field  scattered  from  the  target  is  not  necessarily  the  polarization 
of  the  return  wave  at  the  receive  antenna.  Researchers  such  as  Pilon  and  Headrick 
[24]  have  addressed  this  situation  as  it  applies  to  the  RCS  estimation  problem. 
This  phenomena  must  be  carefully  investigated  when  selecting  the  polarization  of 
the  catalog:  sets. 


CHAPTER  II 


HF  RADAR  SIGNAL  PROCESSING 


2.1  DOPPLER  PROCESSING 


Many  operation  FF  radars  are  used  for  the  detection  and  tracking  of  targets 
on  or  over  ocean  waters.  A  typical  radar  ceil  illuminated  by  a  surface  wave  radar  is 
shown  in  Figure  2.  The  width  of  the  observation  cell  is  determined  by  the  antenna 
bea’nwidth  at  the  radar  operating  frequency,  while  the  length  is  determined  by 
the  time  gating  of  the  received  signal.  To  observe  a  cell  at  a  distance  from  the 
antenna  with  a  cell  length  AR,  the  received  signal  is  gated  in  time  with  values: 


2AR 

T2  =  ri-\ - 


c 


(2.2) 


where  c  is  the  speed  of  light,  and  (t2  —  ti)  is  the  time  gate  used.  With  a  multi¬ 
frequency  radar,  the  beamwidth  of  the  antenna  may  change  since  the  radiation 
pattern  is  a  function  of  freqtiency.  Therefore  the  size  of  the  cell  varies  with  changes 
of  the  radar’s  operating  frequency. 

In  operation,  the  radar  transnuts  a  continuum  of  modulated  pulses  at  some 
given  pulse  repetition  frequency  (PRF).  A  total  of  N  pulses  are  transmitted 
during  the  coherent  integration  time  Tq,  where: 


Therefove  N  coherent  data  records  are  produced  from  the  energy  scattered 
from  the  radar  cell.  These  data  are  Doppler  processed  to  produce  a  power  spectral 
density  estimate  containing  N  data,  with  a  frequency  resolution  and  a  band- 
width  from  to  -f  .  TheoieticHily,  greater  resolution  may  be  obtained 

by  increasing  of  the  coherent  integration  time;  however  changes  in  the  ionospheric 
ray  path,  sea  stale,  and  target  place  a  limit  on  the  practical  integration  time  [22j. 
The  resulting  spectra  shown  in  Figure  3  exhibits  1)  a  return  from  the  target  P'jcif) 
with  induced  doppier  shift  /y,  2)  two  large  clutter  returns  P+3{f)  and  P-sif) 
+fg  and  —fj  respectively,  and  3)  other  relativelj  low  power  clutter  covering  the 
spectrum. 

2.2  DOPPLER  SPECTRAL  ESTIMATE  INFORMATION  CONTENT 

Information  contained  in  the  spectral  estimate  may  be  used  to  estimate  a 
number  of  target  features  for  classification  purposes.  Using  /y,  the  relative  velocity 
of  the  target  may  easily  be  estimated,  giving  the  RTI  designer  knowledge  of  the 
probable  class  of  the  target  (ships,  general  aircraft,  high  speed  aircraft,  etc,).  In 
addition,  estimates  of  the  power  returned  from  the  target  Pxif)  8nd  the  power 
returned  from  some  reference  pRif)  may  be  used  to  estimate  RCS  related  target 
features. 

The  two  large  power  densities  at  -h/jj  and  —fjj  result  from  electromagnetic 
wave  scattering  from  ocean  wave  sets  of  a  certain  wavelength.  The  ocean  waves 
which  excite  this  Bragg  scattering  mechanism  have  a  wavelength  dependent  upon 
the  radar  frequency  [15].  As  explained  in  detail  by  Trisna  [15],  Barrick  [5],  and 
others,  the  amplitude  of  the  Bragg  returns  are  largely  dependent  upon  sea  state. 

Large  variations  in  the  estimates  P^sif)  P-sif)  about  some  means 
(P+Bif))  {P-sif))  ’-Jc  observed  ovei  a  given  observation  period  as  noted 


by  Headrick  and  Pilon  [24].  Therefore  the  Bragg  returns  may  best  be  described 
cis  a  random  process  with  ;ts  mean  being  a  function  of  the  sea  state  at  that  time. 
An  understanding  of  this  process  is  important  if  the  Bragg  lines  are  to  be  used  as 
a  reference  in  the  target  feature  estimation  process. 

2.3  STATISTICS  OF  SPECTRAL  COMPONENTS 

As  discussed,  the  height  of  a  target  return  or  Bragg  line  in  a  given  doppler 
spectra  is  only  an  estimate  of  the  power  returned  from  the  respective  scatter. 
Consider  creating  (M  —  1)  additional  doppler  PSD  estimates  of  the  same  cell  to 
give  a  total  of  M  doppler  PSD  estimates.  Furthermore,  create  these  additional 
estimates  by  observing  the  cell  over  a  continuous  time  span  ioBSi  where 

tOBS  -  ^EST  X  M  (2.4) 

as  shown  in  Figure  4.  Examples  of  real  data  obtained  in  this  manner  may  be  found 
in  Pilon  and  Headrick  [24]. 

Note  in  Figure  4  that  each  spectra  gives  a  different  estimate  for  each  Bragg 
peak  as  well  as  the  target  spectral  peak.  However,  by  averaging  M  samples  taken 
ove.r  M  consecutive  CIT  periods  (k  =  1,...M)  one  may  produce  an  averaged  doppler 
PSD  as  shown  in  Figure  4.  The  averaged  estimated  power  returned  from  the  target, 
approaching  Bragg  wave  set,  and  receding  Bragg  wave  set  are  given  by  Eq.  (2.5), 


Eq.  (2.6),  and  Eq.  (2.7)  respectively. 

M 

(2.5) 

(2.6) 

k=l 
i  M 

P-BUi)=-uT.'’P-BUx) 

k=l 


(2.7) 


If  each  of  the  samples  on  the  right  hand  side  of  Eq.  (2.5)  are  independent  and 
identically  distributed,  the  application  of  the  centred  limit  theorem  as  discussed  in 
Degroot  [10]  gives: 

■  1  M 

liiii  j.  V  =  lim  P^(h)  =  (PtUi))  (2.8) 

ICO  j  Af-*co 

where  {PTifl})  is  the  expected  value  of  the  returned  power  from  the  target  as 
discussed  earlier.  Likewise,  application  on  the  central  limit  theorem  to  Eq.  (2.6) 
and  Eq.  (2.7)  yields  Eq.  (2.9)  and  Eq.  (2.10)  respectively. 

’  M 

T7  E  ‘‘P+Bifl)  =  lim  P^bUi)  =  (P+bUi))  (2.9) 

M— *00  M  M—kx) 

■  1  M 

U  E  '‘P-bUi)  =  lim  P-bUi)  =  (^’-fl(/l))  (2.10) 

Af— *oo  iVi  M—*oo 

Although  this  implies  an  infinite  observation  time  which  is  impractical,  it  also 
implies  that  we  can  continually  improve  the  accuracy  of  our  estimates  by  taking 
more  samples. 

Both  the  CIT  and  the  observation  time  must  be  carefully  chosen  in  order  to 
satisfy  the  requirement  of  independence  and  identically  distributed  samples.  The 
observation  time  must  be  selected  to  cissure  the  identical  distribution  of  samples.  A 
significant  change  in  target  flight  path  during  the  observation  time  may  cause  the 
observed  RCS  to  change  and  therefore  (Pt(/i))  changes  voiding  the  identically 
distributed  criteria.  Likewise  a  change  in  sea  state  during  the  observation  time 
may  alter  {P^Bifl))  {P-Bifl))  same  result.  In  general,  however,  sea 

states  may  remain  unchanged  for  hours  [15].  In  addition,  significant  changes  in  the 
observation  path  (such  as  passage  of  ionospheric  disturbcince  for  sky  wave  radars) 
limits  the  observation  time.  In  fact,  fluctuations  in  the  ionosphere  generally  impose 
the  greatest  restriction  on  the  maximum  observation  time. 


POWER  SPECTRAL  DENSITY 
(ARBITRARY  dB  SCALE) 
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Figure  4;  Averaging  of  iiulcpendcnl  spectral  estimates. 
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There  is  a  limit  to  the  shortest  possible  integration  time  that  produces  inde¬ 
pendent  samples  of  the  Bragg  spectra]  lines.  An  experinientcil  study  by  Baixick 
(13)  concluded  that  an  integration  time  of  25  seconds  or  greater  insures  indepen¬ 
dence  of  spectra,  irregardless  of  sea  state  and  ceil  size,  at  frequencies  in  the  HF 
band  above  7  MHz. 

The  remaining  question  is,  given  a  CIT  of  25  seconds,  what  is  the  relationship 
between  the  standard  deviation  values  of  P^(/i),  and 

As  discussed'by  Barrick  [13],  without  assuming  Gaussian  statistics,  there  is  littie 
to  do  other  than  record  massive  amounts  of  data  and  extract  empirical  results. 
Using  a  Gaussian  assumptions  Barrick  derived: 


_  iia 


(2.1!) 


where  is  the  standard  deviation  of  the  random  variables  pven  in  Eq.  (2.6) 


and  Eq.  (2,7).  Figure  5,  from  Barrick  [13],  shows  the  90%  confidence  interval  lines 
derived  using  the  same  statistical  assumption.  These  curves  ^vc  the  RTI  designer 
insight  into  the  iradeofis  betw^n  the  observation  time  and  the  ^ccu^a^y  of  the 
estimates  P^g{fi)  and  P^(/i).  This  is  important  if  P^g\fi)  or  P^gifi)  are  to 
be  used  in  the  estimation  of  target  feature  for  classification. 


2.4  FREQUENCY  SELECTIVE  FADING 


The  Doppler  spectral  estimate  resulting  from  data  accumulated  via  skywave 
propagation  is  distorted  to  varying  degrees  by  the  Ionospheric  path.  The  distortion, 
referred  to  as  ionospheric  contamination,  is  a  function  of  paran>etcrs  that  may  or 
may  not  be  under  the  control  of  the  radar  operator.  It  is  imperative,  if  data  from 
the  resulting  spectra  are  u?ed  in  an  estimation  and  classification  procedure,  that 


CONFI 


Figure  5:  Relationship  between  deviation  of  average  bragg  line  amplitude  and 
number  of  independent  samples  used  for  averaging. 
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the  distortion  be  rer  /ed  or  at  least  understood  so  it  may  be  compensated  for  in 
the  process. 


Consider  a  target  being  illuminated  by  a  single  wave  which  travels  along  a 
single  ionospheric  path.  The  ionosphere,  which  is  in  constant  motion,  induces  a 
Doppler  shift  on  the  rate  and  direction  of  the  ionospheric  movement.  Now  consider 
illuminating  this  target  N  times  by  waves  which  propagate  by  the  same  layer  whose 
relative  velocity  remained  unchemged  during  the  time  the  data  was  accumulated 


(the  CIT).  The  resulting  spectra  will  be  offset  in  doppler  frequency  by: 

_  2v}  2vj 

hff.„  + 


(2.12) 


where  Vj  is  the  relative  ionospheric  velocity  with  respect  to  direction  of  radar  wave 
propagation  during  the  pass  and  is  the  radar  wavelength.  This  phenomena, 
referred  to  as  Doppler  shifting,  is  easily  detected  in  the  skewing  of  the  spectra  as 
shown  in  Figure  6,  and  can  be  easily  compensated  for. 

Next,  consider  the  case  where  the  CIT  is  relatively  long  such  that  the  relative 
velocity  of  the  propagating  layer  has  undergone  considerable  change.  Eacli  pass¬ 
ing  wave  may  have  a  slightly  different  doppler  component  induced  by  the  quickly 
changing  ionosphere.  The  resulting  contamination  of  the  spectra,  known  a*!  broawi- 
ening,  is  shown  in  Figure  7. 

Spei'tral  broadening  is  highly  tmdesired  for  several  tehsons.  First,  any  reduc¬ 
tion  in  the  amplitude  of  the  target  return  will  m2dce  the  return  more  susceptible  to 
other  noise,  thus  making  any  target  feature  estimated  from  the  data  less  reliable. 
Second,  the  return  from  the  target  may  not  be  significantly  greater  than  the  higher 
order  sea  return  and  noise  surrounding  it,  and  thus  the  tmget  may  be  easily  lost 
in  the  clutter  surrounding  it.  Tliird,  if  the  target  has  a  Doppler  frequency  which  is 


RAL  DENSITY 
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near  a  Bragg  frequency,  th-  resultant  spreading  of  the  Bragg  line  may  essentially 
cover  the  target. 

An  additional  cause  of  Doppler  spreading  is  multipath  propagation  [22].  The 
path  or  paths  the  wave  takes  between  the  transmit  antenna  and  the  target  is  a 
function  of  many  radar  parameters  including  frequency,  launch  angle,  and  polar¬ 
ization,  as  well  as  real-time  ionospheric  conditions.  Different  ionospheric  layers 
may  support  propagation  of  waves  of  the  same  frequency.  This  may  result  in  two 
or  more  waves  with  different  induced  doppler  shifts  both  containing  information 
about  the  target  being  processed  together,  with  a  resulting  spectra  which  exhibits 
spreading. 

A  third  cause  of  ion(»pheric  distortion  restilts  &om  horizontal  spatial  inho- 
mugenity  within  a  given  layer  [22].  A  ray  path  utilizes  an  area  of  a  given  layer 
proportional  in  size  to  the  beamwidih  of  the  antenna.  The  larger  the  beamwidth, 
the  lai^er  the  area  used  for  ref^iion,  and  the  greater  the  effect  of  the  spread¬ 
ing  as  investigated  by  George  and  Marseca  [22].  Since  beamwidth  scales  inversely 
with  frequency,  we  recognize  that  the  effect  of  spatial  inhomogenity  will  decrease 
at  higher  frequencies  for  propagation  via  a  given  layer.  The  BTl  system  designer 
is  interested  in  an  antenna  mth  a  minimal  beamwidth;  however  physical  size  poses 
a  design  restriction.  For  example,  a  beamwidth  of  1°  may  require  antenna  dimen¬ 
sions  on  the  order  of  2  km  [12]. 

In  order  to  use  doppler  spectral  data  in  an  estimation  process,  one  first  desires 
a  measure  of  the  d^ee  of  the  spectral  contamination.  In  general,  the  width  of  the 
Bragg  lines  ^ves  a  reasonable  measure  of  ionospheric  contamination  [22].  A  great 
deal  of  research  has  gone  into  this  area  as  well  as  the  development  of  processes 
designed  to  dean  up  the  distorted  spectra  [22] . 

•An  alternate  approach  is  to  minimize  distortion  by  judicious  selection  of  radar 
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parameters.  The  selection  of  an  observation  time  will  result  in  a  trade-off  between 
resolution  and  broadening.  Since  spectral  peaks  are  largely  uncorrelated  after  25 
seconds  (see  Section  2.3)  we  desire  the  shortest  CIT  above  25  seconds  that  meets 
the  minimum  spectral  resolution  demands  of  the  system.  This  allows  the  greatest 
number  of  independent  spectra  to  be  amassed  in  the  least  total  observation  time. 
Coherent  integration  times  ranging  from  25  seconds  to  100  seconds  are  common. 
The  amount  of  disturbance  a  given  ionospheric  region  is  undergoing  determines  the 
degree  of  spreading  for  a  given  CIT.  Typically,  the  E  and  Eg  regions  are  considered 
most  stable  but  do  not  always  exist.  During  sunrise  and  sunset,  the  ionospheric 
regions  undergo  rapid  change  and  the  degree  of  broadening  significantly  increases. 

The  judicious  choice  of  radar  operating  frequencies  can  also  greatly  reduce 
the  effect  of  ionospheric  contamination.  There  are  however,  important  trade-offs 
between  the  number  of  frequencies  over  which  one  may  obtain  data  and  the  quality 
of  that  data.  The  RTI  system  designer  desires  estimates  of  features  at  a  number 
of  frequencies.  Through  the  use  of  vertical  incidence  ionograms  and  sweep  of 
frequency  backscatter  echoes  one  can  determine  the  frequency  bands  least  likely 
to  result  in  multipath  distortion  [25]. 

One  possible  method  the  RTI  system  designer  may  choose  is  to  use  an  ex¬ 
isting  measure  of  spectral  contamination  with  a  limit,  throwing  away  all  spectra 
contaminated  beyond  that  limit,  and  making  his  estimates  and  classification  de¬ 
cision  based  only  on  the  remaining  data.  This  method  is  introduced  by  Georges 
in  [22].  A  second  approach  is  to  create  feature  estimates  using  all  possible  spectra 
and  then  weight  the  features  in  the  classification  process  according  to  the  degree 
of  distortion  of  the  data  from  which  the  features  were  taken.  A  third  approach 
is  to  estimate  features  from  the  doppler  spectra  whose  values  are  not  effected  by 
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the  spreading.  All  major  peaks  in  the  spectra  are  equally  effected  by  spreading 
such  that  the  ratios  ’(^^01)  •  remain  unchanged  [23]. 

Therefore  use  of  these  ratio  in  feature  estimation  may  provide  estimates  which  are 
relatively  unaffected  by  Doppler  spreading. 


CHAPTER  III 


MULTI-FREQUENCY  PARAMETER  ESTIMATION 

3.1  GENERAL  CONSIDERATIONS 

All  RTI  algorithms  utilize  estimates  of  various  features  of  the  unknown  target 
to  classify  the  target.  Features  such  as  radar  cross  section,  intrinsic  target  phase, 
and  target  speed  may  be  considered  good  taxget  descriptors.  However,  there  is  a 
trade-off  between  the  number  of  features  used  in  a  system,  the  cost  of  the  system 
(processing  time,  etc.),  and  the  resulting  classification  performance.  As  a  result, 
research  has  been  directed  at  finding  the  optimal  types  and  numbers  of  features 
which  result  in  optimal  performance  at  minimal  cost  given  little  or  no  constraints 
in  our  ability  to  obtain  estimates  of  these  features  [2],  [4]  and  [8].  In  an  operational 
HF  radar  cleissification  system  the  numbers  and  types  of  estimated  features  m-%y 
be  severely  limited  due  to  the  constraints  imposed  by  the  KF  channel. 

The  purpose  of  this  section  is  to  study  the  constraints  the  HF  radar  system  im¬ 
poses  on  the  classification  system.  Given  these  constraints,  we  can  then  investigate 
which  and  how  many  of  the  features  available  should  be  used  to  obtain  optimal 
performance.  Previous  research  has  led  to  the  development  of  RTI  algorithms 
which  incorporate  multi-frequency  measurements  of  RCS  and  intrinsic  phase  in 
various  forms  as  classification  features.  The  number  and  selection  of  frequencies 
to  be  used  has  been  studied  with  regard  to  resulting  classification  performance  |2|. 
Different  classification  algorithms  may  place  additional  restrictions  on  frequency 


selection.  However,  in  an  operational  HF  RTl  system,  the  bandwidth  is  likely  be 
constrciined  by  propagation  parameters,  thereby  constraining  the  number  and  type 
of  features  estimated.  The  different  propagation  modes  of  surface  wave  and  sky 
wave  radars  result  in  different  types  of  bandwidth  restraixits  on  the  two  systems, 
therefore  they  are  investigated  independently. 

.3.2  FREQUENCY  BAND  CONSIDERATIONS 

The  frequency  range  of  a  surface  wave  radar  is  not  limited  by  the  ability  of 
the  system  to  propagate  a  wave  over  a  large  range,  but  rather  by  the  propagation 
less  which  iiicreases  with  frequency  [6].  A  target  which  exists  at  some  range  from 
the  radar  site  may  be  illuminated  over  a  frequency  range  of  5  MHz  to  30  MHz; 
however  the  propagation  loss  at  the  upper  frequencies  may  result  in  the  backscatter 
energy  from  the  target  being  on  the  order  of  the  higher-order  sea  clutter  return 
or  galactic  noise  level  therefore  making  detection  and  siny  estimation  of  target 
features  impossible.  In  fact,  the  lower  signal  to  noise  ratio  at  higher  frequencies 
may  result  in  poorer  quality  estimates  of  features  even  though  detection  is  still 
possible. 

Interference  with  other  users  i?i  the  HF  band  may  also  preclude  the  use  of  some 
frequencies.  During  tesnag  of  an  experimental  system  in  Florida,  interference  in 
the  8-13  MHz  band  was  a  siguiticant  factor  in  system  operation  [7].  In  addition, 
backscatter  resulting  from  undesired  £'  -  region  sky  wave  propagation  periodically 
effected  system  performance  at  lower  frequencies  [7].  Therefore  there  are  definite 
limits  to  the  frequency  band  ever  which  estimates  of  features  may  be  obtained. 

Thus  once  system  specifications  have  been  set  (type  of  targets,  minimum 
range,  etc.)  the  frequency  range  for  a  surface  wave  radar  may  be  determined. 
Bands  of  6  MHz  to  8  MHz  may  bt  expected  to  be  used  [7].  From  here  the  RTI 


system  designer  must  determine  which  features  to  estimate,  at  which  frequencies  to 
make  the  estimates,  and  how  many  total  estimates  are  needed  in  order  to  achieve 
maximum  classification  performance. 

The  frequency  band  over  which  an  HF  sky  wave  radar  can  illuminate  a  target  is 
dependent  on  the  current  ionospheric  conditions.  Therefore  recil-time  cissessment 
of  the  ionosphere  is  important  in  order  to  determine  the  band.  Assessment  of 
the  frequency  band  may  be  made  through  the  use  of  backscatter  sounding.  A 
backscatter  sounding  consists  of  radiating  a  frequency-modulated  continuous-wave 
signal  which  is  swept  in  frequency-  then  receiving  the  backscatter  from  the  earth 
after  it  heis  made  two  passes  through  the  ionosphere  [25].  An  oblique  backscatter 
sounding  is  then  produced  by  plotting  the  return  power  level  as  intensity  versus 
time  and  frequency.  An  example  is  given  in  Figure  8. 

A  target  appears  as  a  near  horizontal  line  uu  (he  resulting  plot,  as  shown  in 
Figure  8.  A  direct  measurement  of  the  frequency  band  over  which  the  target  is 
illuminated  may  then  be  made.  Typical  values  range  on  the  order  of  1  MHz  to 
8  MHz. 

It  must  be  noted  that  real-time  monitoring  of  the  ionospheric  path  in  order  to 
determine  the  largest  frequency  band(s)  over  which  a  target  may  be  observed  are 
not  commonplace.  Other  uses  of  sky  wave  radars  such  cis  ocean  waveheight  and 
wind  direction  maps  do  incorporate  real-time  ionospheric  measurements  to  aid  in 
selection  of  frequency  bands  which  are  most  immune  to  multipath  distortion  [21], 
[23],  and  [25]. 

Considerable  research  has  been  and  is  being  conducted  in  the  area  of  real¬ 
time  ionospheric  monitoring.  It  is  apparent  that  these  monitoring  systems  can  be 
of  great  aid  to  an  RTI  system  utilizing  the  HF  sky  wave  radar  as  its  measurement 


3.3  ABSOLUTE  RCS  ESTIMATION 


where  i?  is  the  one  way  path  length  the  wave  travels  to  the  target,  Pj^y  is  the 
average  transmitted  power,  Ct'ja"  is  the  transmitting  antenna  gain,  i^ 

receiving  antenna  gain,  A  is  the  wavelength,  Fp  is  the  one  way  propagation  loss 
factor,  tciT  is  llie  coherent  processing  time,  Ppif)  is  the  received  power  from 
target,  and  £5  is  a  factor  accounting  for  system  losses.  Therefore,  to  directly 
estimate  the  absolute  RCS  of  a  target  at  N  frequencies,  the  above  parameters 
must  be  estimated  at  each  frequenc}-  used.  System  parameters  such  as  P,{v,  G-px^ 
Grx,  tciTt  f  generally  known  or  may  be  accurately  estimated  at 

each  frequency.  The  power  returned  from  the  target  Ppif)  at  each  frequency  may 
be  estimated  from  the  appropriate  Doppler  spectra  as  discussed  earlier.  However, 
the  range  {R)  and  propagation  loss  (Fp)  may  be  much  more  difficult  to  estimate 
depending  on  the  propagation  mode  involved. 

For  a  surface  wave  radar,  the  range  to  the  target  may  easily  be  estimated 
through  time  gating  or  other  signal  processing  strategies.  The  propagation  loss, 
however,  is  a  function  of  the  sea  state  over  the  path  as  well  as  frequency  and 
thus  becomes  more  difficult  to  accurately  estimate;  see  Norton  [6j.  A  sky  wave 
radar  poses  a  problem  when  attempting  to  estimate  the  path  length  R.  The  first 
difficulty  results  from  path  length  being  a  function  of  the  ionospheric  structure 
which  is  constantly  changing  with  time.  In  addition,  waves  of  difl'eiejit  frequencies 


may  travel  different  paths  to  illuminate  the  same  target  thus  making  R  not  only 
a  function  of  time  but  also  of  frequency  [14j. 

To  circujuvent  the  difficulties  of  estimating  Fp  in  a  surface  wave  radar  and  R 
in  a  skywave  radar,  a  reference  may  be  employed.  Possible  references  include  land 
masses,  ocean  markers,  or  Bragg  line's.  Knowledge  of  the  radar  c  'oss  section  of  the 
reference  then  allows  estimation  of  the  target  RCS  at  each  frequency.  .4  method 
of  calibrating  Bragg  lines  for  use  as  an  absolute  reference  is  given  in  Trisna  (1.5j. 

3.4  RELATIVP  RCS  ESTIMATION 

Bragg  lines  may  also  be  employed  as  a  relative  reference.  Figure  9  shows  two 
spectra  obteiined  at  two  frequencies  /j  and  f-y-  Knowing  the  expected  change  in 
the  mean  of  the  Bragg  line  amplitude  with  frequency  (which  is  cissumed  to  be  zero 
in  this  example)  the  relative  RCS  estimates  of  the  target  can  be  constructed  as 
shown  in  Eq.  (3.2)  and  Eq.  (3.3). 

(rjifi )  =  (20  -  Z0)dB  =  -lOdB  (3.2) 

=  (15  -  Z0)dB  =  -20dB  (3.3) 

Therefore,  the  RCS  of  the  target  at  fi  is  1/lOth  the  RCS  of  the  reference,  and  the 
RCS  of  the  target  at  /2  is  1/lOOth  the  RCS  of  the  reference.  Thus  the  relative 
change  of  the  target  RCS  has  been  estimated  without  absolute  calibration  of  the 
Bragg  line.  Therefore,  the  use  of  the  relative  RCS  of  the  target  as  a  possible  feature 
for  classification  purposes  should  be  evaluated. 

3.5  RELATIVE  PHASE  ESTIMATION 

Previous  KTI  studies  have  shown  thal  intrinsic  taiget  phase  often  proves  to 
be  a  powerful  feature  when  used  for  clHs.^ificatio;i  [21,  (81.  However,  difficulties  in 


.  DEHSITY  POWER  S?ECTRAL  DENSITY 

SCALE)  (ARBITRARY  d3  SCALE) 


^ccuralcly  esiimating  the  path  ienglh  Jo  the  target  a,s  in  flu  <'a.sc*  of  thf*  sky  wave 
prcj>agatiun  nia.y  preclude  the  accurate  e$lima(iou  of  intrinsic  '•'rget  phase.  There¬ 
fore  an  alternate  niethof!  of  utilizing  target  phase  inforniatioi*  has  been  ffudied. 

('o)!si(b»r  the  pliases  of  two  radar  signals  at  fre(juency  /j  an«';  /2.  whir)-,  ar." 
returned  from  the  same  tp.rget  as  gix'en  by  Eq.  (3.4)  and  Eq.  (3.5|. 
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where  i.,  the  round  trip  path  length  of  wave  at  frequency  /;,  R j.,  is  the  round 
trip  path  length  of  wave  at  frequ- ncy  ^2  —  T^i  ^/,  ‘"•^1 

the  intrinsic  taigef  pliases  at  /j  and  /•>  respectively. 

By  multiplying  Eq.  (3.4)  and  Eq.  (3-0)  by  A]  and  A-j  respectively  and  thcr 
.suijtracting  E<|.  (3.4)  from  Eq.  (3.-5),  vve  are  left  with: 


(3.0) 


where  [Rf,  -  Rf^}  is  the  difference  between  round  trip  pa»h  lengths  at  /2  *=’-'** 
f\  respectively.  Through  careful  .selection  of  A/  -  J2  ~  fi  »‘«>d  using  if>»iogran-s 
to  estimate  the  present  structuie  of  Ihe  iono:,pi;ere,  it  jiiay  be  pos.s:b|e  1:>  accu- 
raicly  estimate  {Rf^  Rf^ ).  By  multiplying  this  estimate  hv  2a  radians,  -an-l  t},eu 
subtracting  this  result  from  Eq.  (-3.t3).  we  are  tefl  with  Eq.  (3.7). 


(.3.7) 


Further  division  <>v  givc-s  Eq.  !3.-'':.  winch  for 
appro.ximaled  l>y  Kq.  {3.*1). 
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Therefore  If  is  a  feature  dependent  on  the  change  of  the  characteristic  piiase  of 
the  target  associated  with  fo  and  /]. 

Figure  10  demonstrates  a  physical  interpretation  of  the  IT  function  for  a 


Note  ihe  If  function  is  a  measure  of  the  distance  between  the  scattering  for  the 
two  frequencies,  and  not  the  absolute  po.sifion  of  the  scattering  point  at  any  one 
frequency. 
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A  niults-freqxiency  radar  syste*«  iihiminating  a  target  with  M  total  adjacent 
frequencies  could  generate  (M  —  1)  If'  function  estimates.  These  estimates  could 
then  be  employed  as  additional  features  to  help  identify  the  target.  In  general,  the 
larger  the  A/,  the  more  the  patiis  taken  bj'  the  t.wo  waves  differ,  and  the  harder 
it  is  *0  a<-curately  estimate  (Hj,  —  Rf^)-  Ho\'-ever,  large  A/  values  would  tend  to 
excite  scattering  points  on  the  ziircraft  which  are  more  separated  in  distance,  pos¬ 
sibly  leading  to  IT  functions  which  contain  more  information  about  the  individual 
targets.  Therefore,  our  ability  to  accurately  estimate  the  change  in  path  lengths 
with  frequenct'  for  a  given  ionospheric  condition  may  determine  our  mzoumum  A/. 
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Cl  [AFTER  IV 

RADAR  SYSTEM  SIIV  ULATION  AND  EVALUATION 

STUDY 

4.1  INTRODUCTION 


t4 


In  order  to  study  the  perforr.  ance  of  possible  radar  target  classification  sys¬ 
tems  wliich  utilize  HF  radars  as  a  measurement  device,  a  simulated  system  has 
been  developed.  Over  the  past  several  years  the  ElectroScience  Ivaboralory  has  de¬ 
veloped  a  large  multi-freciuency  data  base  consisting  of  aircraft,  ship,  and  ground 
vehicle  radar  signatures.  This  data  base  has  been  u.sed  to  investigate  several  RTI 
methods  for  the  various  classes  of  targets.  A  data  base  con-sisting  of  calibrated 
monostatic  radar  returns  of  •'5  aircraft  was  chosen  for  the  experimental  portion  of 
this  study.  Although  classification  of  the  ships  and  ground  vehicles  ij  of  consid¬ 
erable  interest,  the  large  amount  of  j)revious  work  with  the  aircraft  data  base  [2] 
provides  a  strong  foundation  of  previous  work  for  comparing  results. 

The  ElectroScience  Laboratory  compact  radar  range  facility  was  used  to  cre¬ 
ate  the  data  base.  A  basic  block  diagram  of  the  measurement  system  and  a  picture 
of  the  actual  facility  is  given  in  Figures  11  and  12  respectively.  A  more  detailed 
explanation  of  the  .system  may  be  found  in  [3j.  In  general,  the  meastireiiienl  facility 
simulates  the  illumination  of  a  target  by  a  radar  wave  which  is  approximately  plane 
in  the  region  of  the  target  (locally  plane).  This  would  be  true  of  any  HF  radar 
system.  The  complex  backscattcred  fields  from  tlie  target  and  other  undesired 
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scatters  are  then  measured  over  a  frequency  band  of  1-12  GHz.  To  obtain  the 
high-accuracy  range-independent  complex  signature  of  the  target,  a  calibration 
sequence  described  in  detail  in  Kimball  (I7j  was  employed.  Further  details  con¬ 
cerning  data  scaling,  windov.'ing,  and  other  signal  processing  details  are  given  in 
Harris  [26],  and  Kimball  [17j.  Five  metallic  scale  models  of  the  Concorde,  DCIO, 
707,  727,  and  747  were  used  as  targets.  Silhouettes  and  dimensions  of  the  real 
aircraft  are  given  in  Figures  13  through  17.  Specifications  of  the  -5  aircraft  data 
base  are  given  in  Table  1. 

4.2  RADAR  SYSTEM  SIMULATION 


A  flowchart  of  the  Radar  System  Simulation  and  Evaluation  (RSSE)  program 
is  shown  in  Figure  18.  The  program  begins  with  selection  of  N  complex  radar 
signatures  of  the  5  targets  at  frequencies  •  m/jv?  where  the  frequencies  are 

chosen  by  the  operator  to  represent  those  which  may  be  encountered  in  a  typical 
HF  radar  system.  In  addition,  one  azimuth  angle  and  one  polarization  is  selected, 
thus  assuming  these  parameters  have  been  accurately  estimated.  The  target  data 
base  is  then  accv^ssed,  the  data  are  appropx'iately  scaled  [8],  and  two  sets  of  data 
denoted  the  test  and  catalog  set  are  created.  Each  element  of  the  two  sets  may  be 
represented  by: 

+  ,A'i.(/i)  (4.1) 

where 


Rk{fi)  =  Ak(fi)cos(e^.{fi)),  (4.2) 

A');.(/,)--4i.(/,-)sin(0j^.(/,)).  (4.3) 


and 

33 

I 


I 


I  .<71  tfT  yy  J,  j:XjrXi<TUni.>»  Jj*'^jC%>i.-«Ma'VU^. 


rfcF.rjvir-ijC-'. 


Figure  11:  Compact  range  block  diagram,  (from  Kamis  [2j), 
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Figure  12:  OSU  ESL  compact  range,  (frojn  Kamis  j2i). 
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ire  17:  747  silhouette  (from  Kamis  [2]). 
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Low-IVequency  Data  Base 

Frequency  formatted  data  strings  from  1  to  12  GHz 
Availible  polarizations: 

Transmit  Horizontal,  Receive  Horizontal  (HH) 
Transmit  Vertical,  Receive  Vertical  (VV) 

TVansmit  Horizontal,  Receive  Vertical  (VH) 

Availible  Azimuth  angles  6  Elevation  =  0°: 

0®  to  180°  by  10°  end  15°  increments 
Common  aircraft  bandwidth: 

7.6  -  60  MHz  “Scaled” 

Maximum  number  of  usable  frequencies:  209 


Table  1:  Data  ba.se  speeificatious  (from  Kaniis  [2]). 
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■^kifi)  -  square  root  of  the  radar  cross  section  of  target  at  frequencj' 


A  corrupted  ffature  vector  representing  an  unknown  target  was  constructed 


by  first  generating  2N  total  noise  samples  taken  from  an  independent  Gaussian 


zero-mean  random  process  with  noise  power  P;v'  ^’^cn  by  : 


PjV  =  101ogjo(  — IdPsrn 


where  a  is  the  standard  deviation  of  the  process  in  meters  (3J.  The  iioise  samples 


were  then  added  to  the  N  high  accuracy  complex  radar  signatures  of  a  known 


target  in  the  test  set  to  create  A’  corrupted  complex  data  of  the  form: 


Ru(fi)  =  RT{fi)  +  N{fi) 


A'«(/i)  =  AT(/,)-f  A''{//) 


w'here  Pj-I/, )  is  defined  in  Eq.  (4.2)  for  target  T  in  the  test  set,  A'^l/i)  is  defined 
in  Eq.  (4..3)  for  target  T  in  the  test  set.  and  A'(/,),  N'(fi)  are  samples  from  the 


independent  zero-mean  gaussian  noise  process. 


The  resulting  A*  complex  corrupted  data  which  may  be  represented  by 


RAfi)  rjXAfi)  =  A(fi)ej^"^Rl  i  =  h...,N 


constitute  the  unknown  target  feature  vector.  This  corrupted  feature  vector  is 


then  tested  agciinst  the  feature  vectors  of  the  targets  in  the  catalog  set  by  means 


of  a  classifier  function  as  shown.  A  Euclidian  distance  metric  algorithm  commonly 


referred  to  as  a  Nearest  Neighbor  (NN)  algorithm  was  emi)loyed  the  classi¬ 


fier  function.  The  NN  algorithm  computes  a  distance  |d„j^.)  as  a  measure  of  the 


distance  between  the  unknown  target  »  and  the  catalog  target  k\  as  shown  in 


Eq.  (4.3). 


.".<•'1  tc*.*.*-.-*.- -csc.-tt  »»“K?«*aac*irr.  *4^-i'''t‘'rscTiarsc;r' 
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wh'^rp  !  is  the  feature  of  larg^’t  k  at  frequency  f,;  extracted  fiom  t!i?  catalog 
set.  and  Tuif,)  is  tlse  featuie  of  unknov.-ii  target  ii  at  frequency  /,  extracted  from 
the  corrupted  fea.iure  vector. 

After  tcstine  against  the  appropriate  features  of  all  5  possible  catalog  targets 
(k  -  I  to  5)  the  unknown  target  is  classified  as  the  catalog  target  which  produced 
the  smallest  distance  measure.  Whether  the  classification  attempt  was  succesful 
or  not  IS  known  since  the  identity  of  the  unknown  target  is  known  a-priori.  This 
simulation  process  is  lepeated  100  tin.c.s  for  each  test  target,  with  independent 
noise  being  added  each  time  to  create  the  simulated  feature  vector.  Results  of  a 
run  for  a  given  noise  power  are  plotted  as  a  confusion  matrix  as  shown  in  Figure  19. 

TT.e  average  uuscleissification  percentage  has  a  margin  of  error  since  it  is  the 
result  of  a  finite,  and  not  an  infinite,  number  of  independent  experiments.  Details 
of  the  error  margin  may  be  found  in  Degroot  jlOj.  In  general,  the  error  margin  is 
on  the  order  of  2  to  -3  percent  for  -5  targets  and  iOO  experiments  |2j. 


4.3  CHANNEL  SELECTION  AND  SIMUL.ATION 


in  general,  tiie  ability  to  accurately  estimate  target  features  such  as  abso¬ 
lute  KCS  and  intrinsic  ;ohase  may  be  described  by  defining  channels  with  various 
rr.aractcristics  wdiich  affect  the  complex  barkscattcccd  coelncicnts  of  the  targets. 
This  section  describes  four  .'-hannels  wish  dilferent  characteri.-'tics  which  allow  dif¬ 
ferent  targes  features  to  be  estimated.  By  appropriately  proccs.siii.c,  Uu-  complex 
backseat tered  components  from  the  test  set  foe  each  chnuuel  model,  and  ah'T- 
ins  the  complex  bac.kscatl-’r  component.*:  of  the  catalog  set  if  appropriate,  the 
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Figure  19:  Confusion  matrix  (from  Kaiuis  |2]). 
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resulting  classificalioii  performance  given  the  coustreiir.t?  of  each  channel  may  be 
synthesized. 

4.3.1  Coherent  Gaussian  Channel  Model  Simulation 

A  channel  which  adds  coherent  independent  zero  mean  Gaussian  noise  to  the 
complex  backscatter  coefficients  of  target  T  is  shown  in  Figure  2C.  Note  tha.t 
both  the  amplitude  and  phase  of  the  high-accuracy  backscatter  terms  ha-ve  been 
corrupted.  The  resullirig  estimates  Rail,)  Xuif,)  lepresent  th«*  corrupted  real 
imaginary  parts  of  the  complex  backscatter  coefficients  of  an  unknown  target  to 
be  classified. 

The  impleiiientation  of  the  Channel  /  model  into  the  RSSH  program  ha-s 
previously  been  described  in  Section  4.2.  .An  algorithm  denoted  the  Coherent 
Neaj  st  Neighbor  .Algorithm  given  by: 

r,v  -ji 

'iu.L-  -  T,^Rh{fi)  -  ««{/.))-  +  (.Vtt/i)  -  .Y„(/,))-  (4.9) 

.*=1 

where  Ri{f,),  .Yjtl/i)  are  the  complex  backscatter  coefficient  of  catalog  target  k. 
and  Ruift )-  -Yu(/j )  are  the  simulated  real  and  imaginary  parts  of  complex  backscat- 
ter  coefficient  of  unknowzt  target  u,  is  then  used  to  classify  the  unknown  target  V 
as  being  the  catalog  target  A*'  such  that  d„  |.»  is  minimized.  It  is  important  to 
note  that  intrinsic  target  phase,  although  not  directly  estimated  at  the  output  of 
Channel  1.  is  utilized  in  the  estimation  process. 

4.3.2  Non-Coherent  Gaussian  Channel  Model  Simulation 

.A  channel  which  corrupts  the  an:p!studcs  of  the  complex  backscatter  coeffi¬ 
cient. s  throuuh  coherent  addition  of  independent  zero  mean  Gaiussian  noi>e.  while 
removiiie  all  intrinsic  phase  information  is  shown  in  Figure  21 .  This  channel  models 
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a  radar  system  where  accurate  esthiiation  of  intrinsic  target  phase  is  not  feasiijle. 
As  a  result,  only  estimates  cf  backscatter  amplitudes  are  available  for  classilication. 

Tlse  implementation  of  the  Channil  2  model  into  the  RSSE  program  has 
previously  been  described  in  Section  4.2.  .An  algorithm  denoted  the  Nou-t  onerent 
Nearest  Neighbor  .Algorithm  given  by: 


<^u,k  - 


r  A' 


(4.10) 


•.vhere  Ai-if,)  is  the  high  accuracy  amplitude  of  complex  backscatter  c«>efficient  of 
catalog  target  k\  and  Ru(fi)  is  the  simulated  amplitude  of  complex  backscatter 
coefficient  of  unknown  target  «,  is  then  used  to  classify  the  unknown  target  a  as 
being  the  catalog  target  k*  such  that  i-*  is  minimized. 


4.3.3  Multiplicative  Component  Channel  Model  Simulation 


.A  channel  which  ?iiuUiplies  the  amplitudes  of  an  unknown  target's  backscatter 
coefficients  by  an  unknown  constant  M  in  addition  to  adding  noise  and  renjoving 
phase  information  is  shown  in  Figure  22.  This  channel  models  a  radar  system 
where  accura.te  estimation  of  absolute  RCS  is  not  practical  due  to  the  lack  of  a 
calibrated  reference,  but  estimation  of  relative  RCS  as  described  in  Section  -3.4  is 
feasible. 

One  possible  method  of  using  -4^^(/,-)  data  out  of  Ckanncl  S  involves  the 
removal  of  ihe  unknown  factor  M.  To  achieve  this,  the  average  power  of  the 
corrupt  vector  cut  of  Channel  7  is  calculated  by: 
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Figure  20:  Coherent  AWGN  channel  model. 
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Figure  21:  Non--Cohcrent  AWGN  channel  model. 
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=  M-P 


(4.13) 


where  P  is  the  average  power  of  the  noise  vector  if  M  —  }.  A.  new  corrupt  vector 


with  values  given  by: 


AUfi) 


s/W 


(4.14) 


.<(.f,) 


sliRrUi)  +  NUi)9^{XT{fi)  +  N>[fi)f^ 


(4.15) 


may  then  be  constructed.  This  resulting  vector  now  has  values  independent  of  M 
while  maintaining  tii  ^tio’s  of  a.nplitude  estimates  at  different  frequencie.^.  i.e.: 

^  for  1  <  i  <  JV,  1  <  )  <  N.  (4.16) 

The  A^(fi)  features  may  now  be  used  for  target  classification  with  a  properly 


1  <  .?  <  N- 


(4.1fi) 


chosen  catalog  set. 


The  simulation  of  utilizing  features  passed  through  a  channel  with  an  unknown 
multiplicative  component  first  required  normalization  of  the  catalog  set.  The  power 
(or  average  RCS)  of  catalog  target  k  is  given  by  Eq.  (4.17). 


h  =  i 
^  1=1 


(4.17) 


Dividing  each  data  of  catalog  target  k  by  the  square  root  of  Pf.  gives  a  new  set  of 
data,  where  the  magnitude  of  each  point  in  the  vector  is  given  by  Eq.  (4.18). 


AR(f  Af.(fi) 


(4.18) 


The  new  data  A^(f,)  represent  the  relative  radar  cross  section  features  of  catalog 
target  k  which  will  be  used  in  the  classification  routine.  Note  that  the  latio  of  radar 
cross  section  of  any  two  frequencies  /,  and  fj  for  catalog  target  k  is  pre-served  as 


shown  in  Eq.  (4.19). 
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Figure  22:  Multiplirative  aWGN  channel  model. 
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1  <  j  <  iV  ,  1  <  7  < 


(A-(/,)P  _  (/if(/i))- 
(^i.(/,))“  (>tf(/i))- 

The  normalization  described  in  Eq.  (4.17)  and  Eq.  (4.18)  is  carried  out  on  ali 
catalog  targets  (k  =1  to  .5).  As  a  result,  each  catalog  target  has  an  average  power 
(or  average  RCS)  given  by: 
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Therefore,  the  overall  average  size  of  each  target  in  the  catalog  has  been  removed. 

The  normalization  is  graphically  depicted  in  Figure  23.  The  complex  backscat- 
tcred  amplitudes  at  frequencies  /j  and  /2  for  two  catalog  targets  are  depicted  by 
vectors  Vj  and  V2*  Vector  V2  represents  a  relatively  large  target  with  an  average 
RCS  on  the  order  of  1000  sm,  and  V'j  represents  a  smaller  target  with  an  average 
RCS  on  the  order  of  400  sm.  Note  the  distance  between  the  two  points  in  space 
representing  the  targets  is  on  the  order  of  10  to  15  meters.  After  normalization, 
the  two  targets  are  represented  by  vectors  Vi  and  Vo-  Note  the  distance  between 
the  points  in  space  representing  the  two  targets  has  significantly  decreased. 

The  addition  of  coherent  i.i.d.  noise  to  the  complex  backscatfered  components 
of  a  te.st  target  to  form  a  corrupt  test  vector  was  carried  out  as  described  earlie». 
After  the  noise  was  added  the  power  of  the  test  vector  wa.s  calculated  by: 
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producing  the  normalized  amplitudes  of  the  backscatter  estimates  of  the  unknown 
target.  This  process  would  remove  any  unknown  multiplicative  component  from 
the  test  vector  as  shown  earlier.  It  is  critical  that  the  noise  be  added  before  the 


normalization  process  is  completed  in  order  to  directly  compare  the  relative  RCS 
to  the  absolute  RCS  as  possible  claorificaiiou  features.  If  normalization  is  done 
before  noise  is  edded,  the  signal  to  noise  ratio  will  be  much  lower  resulting  in 
significantly  poorer  classification  performance.  The  unknown  target  may  then  be 
tested  against  each  of  the  catalog  targets  by  use  of  a  Nearest  Neighbor  Algorithm 
of  the  form: 

'  N  1  ^ 

=  E  '  (4.26) 

.1=1 

where  and  are  given  in  Eq.  (4.18)  and  Eq.  (4.25)  respectively,  and 

the  identity  of  the  unknown  target  is  chosen  as  the  catalog  target  k*  for  which  d,j  jt 
is  minimized. 

4.3.4  Additive  Component  Channel  Model  Simulation 

Since  the  case  of  using  a  channel  model  with  an  unknown  multiplicative  factor 
M  is  of  interest,  it  is  natural  to  extend  the  problem  to  the  ca.se  of  an  unknown 
additive  component  C  as  shown  in  Figure  2-i.  The  addition  of  the  unknown  com¬ 
ponent  C~  under  the  radical  sign  as  shown  represents  the  estimation  of  RCS  of  a 
target  at  2  oi  more  frequencies,  where  only  the  diflcrence  in  RCS  is  known.  Ff)r 
example,  the  difi'cience  in  RCS  of  a  target  at  f\  and  /2  may  be  known  to  be  50G  sm 


but  whether  the  DCS  at  fi  and  fo  is  10  sm  and  510  sni  respectively,  or  1000  sm 
and  1500  sni  is  not  known.  Although  this  problem  may  or  may  not  arise,  it  is  of 
considerable  interest  to  the  RTI  designer  to  determine  the  achievable  classification 
performance  given  a  channel  with  this  constraint. 

To  remove  the  unknown  additive  component  C'“,  the  average  j)ower  of  the 
corrupt  vector  at  the  output  of  Channel  4  is  calculated  by: 


(4.27) 


P^'  =  C-  +  P 


(4.28) 


where  P  i.s  the  average  power  if  C"  =  0.  The  vector  is  then  normalized  Ijy: 
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A? (/.)  =  [(ft; (.f.)  +  +  (Xrifi)  d  NU))' 

d  O'”  [Pnorm  r(r--P)j]*  (4.30) 

A{!if,)  =  [(ftr(/  )  d  iV(/.))-  +  (A*r{/.)  +  A’'(/.))2  +  [P  -  P, (4.31) 

giving  a  new  vector  independent  of  6’,  while  maintaining  the  difference  in  power 
(RCS)  between  data  at  any  two  frequencies. 

Features  constructed  from  a  channel  with  an  unknown  constant  additive  com¬ 
ponent  are  referred  to  as  relative  difference  RCS  (RDRCS)  features.  The  features 
were  implemented  into  the  RSSE  program  by  reco:istructing  each  target  in  the 
catalog  set  so  its  new  average  power  is  the  same  for  all  targets  while  the  dif¬ 
ference  between  RC^S  at  a  given  fiequencv  and  the  average  RCS  is  maintained.  The 
same  reconstruction  of  (he  corrupted  test  vector  remove.s  ai  y  unknown  additive 
component  from  the  e.stimatc.s.  .A.s  .i  resuit.  each  target  in  tlie  catalog  .set  is  given 
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Figure  24:  Additive  AWGN  channel  model. 
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the  siinip  avf'rage  RCS,  with  the  difference  between  the  RCS  at  each  frequency  and 
the  average  RCS  being  maintained. 

The  normalization  of  the  catalog  set  in  the  RSSE  program  was  accomplished 
as  de.scribed  above.  Since  the  amplitude  of  the  complex  backscatter  coefficients 
in  the  catalog  set  is  actually  the  square  root  of  the  radar  cross  section  at  each 
frequency,  tlie  amplitudes  were  reconstructed  using; 
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=  [(-4a.(/;))“  +  (P„or„.  -  Pji.)]' 


(•1.32) 
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An  example  of  the  normalization  effect  on  two  catalog  targets  is  given  in 
Figure  2-5.  Consider  target  1  in  the  catalog  set  with  RCS  values  SOOsm,  400sm,  and 
300sin,  and  target  2  with  RCS  values  of  GOOsm,  lOOOsm,  and  1400sm  at  frequencies 
/l,  f-y  and  /;;  as  shown.  The  average  RCS  over  the  three  frequencies  is  400sm  and 
lOOOsm  for  targets  1  and  2  respectively.  After  the  normalization,  the  average  RCS 
of  both  targets  is  IGOOsm  (Puorm  =  IGOOsm)  as  shown.  As  a  result,  absolute 
target  size  is  lost;  however,  the  difference  between  the  RCS  at  any  frequency  and 
the  average  RCS  is  maintained. 

Coherent  noise  was  added  as  always  to  the  complex  backscatter  coefficients  of 
target  T  of  the  test  set  to  create  the  noisy  feature  vector.  The  average  pov/er  of 

the  noisy  feature  vector  wcis  then  calculated  as: 
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The  amplitudes  of  the  noisy  feature  vector  were  then  reconstructed  as: 
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Figure  25:  Normalization  of  two  targets. 
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These  data  represent  the  estimates  of  the  aniplititdes  of  the  complex  coetlicients 
passed  through  channel  4  after  the  additive  unknown  component  has  been  removed. 


As  a  result,  the  Nearest  Neighbor  Algorithm  given  by: 


where  1  and  are  given  by  Eq.  (4.32)  and  Eq.  (4.35)  respectively,  wa.s 

used  to  c;a.ssify  the  unktiowa  target  »/  as  the  catalog  target  ic*  for  the  b*  such  ♦hat 
minimized. 

It  is  important  to  note  that  although  the  power  of  the  catalog  vectors  and 
c<irrupted  vector  is  altered,  the  signal-io-noise  ratio  of  the  initial  corrupted  vector  is 
not  chaiiged.  and  thus  classification  perfoiinance  will  not  be  function  of  the  power 
level  chose:i  for  normalization.  Simulations  were  repeated  using  normalization 
power  levels  of  4000  sm,  6000  sm.  and  8000  sm  to  insure  this  result. 

4.4  ASSESSMENT  OF  MISCLASSIPICATION  CURVES 

A  t\  pical  set  of  misclassificafion  curves  is  shown  in  Figure  26.  Performance  of 
a  given  curve  with  respect  to  another  may  be  related  by  the  relative  noise  power 
immunijy.  The  noise  power  immunity  of  a  curve  is  determined  by  tiie  difference 
between  the  noi-se  power  of  the  curve  and  the  reference  curve  evaluated  at  a  con¬ 
stant  misclassifiration  level  (generally  between  5  and  40  percent ).  The  lu  ise  power 
immunity  is  an  approximate  measure  of  the  additional  noise  power  which  must  cor¬ 
rupt  tin-  according  data  to  produce  a  curve  equivalent  to  the  reference.  A  positive 
noise  power  immunity  indicates  better  classification  performance  with  resperl  to 
the  reference,  wliile  negative  noise  power  immuiiity  indicate.^  the  opposite. 


CHAPTER  V 


RSSE  RESULTS  FOR  AWGN  CHANNELS 

5.1  Frequency  Band  Study 

The  frequency  band  over  which  a  target  at  a  given  range  may  be  iiluminaf  i.-d 
is  limited.  As  a  result,  only  estimated  features  from  the  limited  band  are  available 
for  classification.  Thus  a  study  has  been  conducted  to  find  the  effect  ot  frequency 
band  size  on  classification  performance  for  the  four  channels  described  in  Sections 

4.3.1  througN  4.3.4. 

Figure  26  shows  the  change  in  performance  as  the  frequency  band  of  channel  1 
is  decreased.  Note  the  overall  performance  is  not  significantly  effected.  Figure  27 
demonstrates  the  greater  loss  of  performance  with  decreasing  frequency  baud  when 
channel  2  used.  A  significant  loss  of  noise  power  immuni'y  of  approximately  8  dB- 
sm  is  encountered  when  the  frequency  band  decreased  from  4  MHz  to  2  MHz. 

Similarly,  Figure  28  demonstrates  the  even  greater  decrezise  in  classification 
performance  encountered  by  decreasing  the  frequency  band  when  channel  3  int  re¬ 
duces  an  unknown  multiplicative  component.  Losses  on  the  order  of  '5  md  12  dB- 
sm  noise  power  for  a  given  misclassificafion  percentage  appear  as  the  frequency 
band  is  decreased  from  8  MHz  to  4  MHz,  and  then  4  MHz  lu  2  MHz.  F'gitre  2‘'i. 
which  was  generated  utilizing  the  unknown  additive  component  channel,  expvess<  s 
the  most  significant  decrease  in  perfo’  inance  as  the  band  i.‘i  decreased  from  S  Mliz 
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Figures  26  through  29  demonstrate  some  qualities  unique  to  ea~h  of  the  four 
channels.  Channel  1  demonstrated  a  very  small  change  in  performance  compared 
witii  the  three  other  rhannels  as  the  frequency  band  was  signincantly  decreased. 
Botli  channel  2  and  channel  3  showed  a  much  larger  degradation  in  performan''e 
when  the  band  was  reduced  to  2  MHz  as  compared  to  the  reduction  to  4  MHz. 
This  impliec-  a  frequency  band  of  not  less  than  4  MHz  would  be  highly  desirable 
g»\en  a  '■hannel  which  does  not  allow  the  estimation  of  phase  information.  This 
also  implies  that  a  doubling  of  the  frequency  band  from  4  MHz  to  S  MHz  does  not 
provide  the  same  degree  of  improvement  in  classification  perfornrance.  Therefore 
the  cost  of  having  the  ability  to  obtain  features  over  an  8  MHz  frequency'  band,  such 
as  a  second  radar  site,  may  not  be  warranted.  Conversely,  channel  4  demonstrated 
significant  performance  degradation  in  both  czises  of  band  reduction.  Therefore, 
given  a  channel  with  an  unknown  additive  component,  a  significant  increase  in 
performance  may  be  obtained  by  extending  the  frequency  band  beyond  4  MHz. 
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Parameter  Putlcr  Observation:  Frequency  Banil  Size 

CMiannei  Modeb  Coherent  AWGN  Channel 

Azimuth  angle:  45 

Start  Frequency:  8-0  KIHz 

Slop  Fre<|uency:  Varies 

Nuir.ber  of  Freesuennes:  10 


Ng6E  0)6-31) 

8-0  MHz  to  10.0  MHz  Band 
8.0  MHz  to  1?  0  MHz  Band 
8.0  MHz  to  13.0  MHz  Band 


Figure  25:  C^oherent  channel  band  size  study. 
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Parameter  Under  Observation:  Frequency  Band  Size 

Channel  Model;  Multiplicative  Component  AWGN  Channel 

Azimuth  angle:  45 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  Varies 

Number  of  Frequencies:  10 


_ 8.0  MHz  to  10.0  MHz  Band 

- 8.0  MHz  to  12.0  MHz  Band 

—  -  —  8.0  MHz  to  16.0  MHz  Band 


F'igiuv"  28:  Multiplicative  component  channel  band  size  study. 
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Parameter  Under  Observation:  Frequency  Band  Size 

Channel  Model:  Additive  Component  AWGN  Channel 

Azimuth  angle:  45 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  Varies 

Number  of  Frequencies:  10 
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8.0  MHz  to  10.0  MHz  Band 
8.0  MHz  to  12.0  MIJz  Band 
8.0  MHz  to  16.0  MHz  Band 


Figure  29:  Additive  component  channel  baud  size  study. 


5.2  Number  of  Features  Study 


A  study  was  conducted  to  determine  the  relationship  between  the  number  of 
features  used  and  file  classification  performance  for  the  four  channels.  Previous 
research  [2]  concluded  that  the  minimum  frequency  sampling  interval  which  will 
satisfy  Shannon’s  sampling  criteria  is  given  by: 


(5.1) 


where  c  is  the  speed  of  light  and  L  is  the  length  of  the  target.  Therefore  given  an 
average  target  length  of  50  meteis,  one  finds  approximately  2,  3,  and  5  frequency 
samples  a:e  required  to  me^t  the  above  criteria  for  band  widths  of  2  MHz,  4  MIIz, 
and  3  MHz  respectively.  These  values  were  used  as  minimum  guidelines  to  be  met 
for  each  corresponding  frequency  baud.  Since  Kamis  [2]  examined  the  effect  of  the 
number  oi  features  at  large  frequency  bands  >  8  MHz)  for  channels  1  and 

2,  these  bands  for  these  channels  were  not  included  in  the  study. 

Figme  30  shows  a  significant  increase  in  noise  immunity  of  approximately 
5  dB-sm  as  the  number  of  features  through  channel  1  are  increased  from  3  to  .5 
to  10  over  the  2  MHz  band.  The  increase  in  noise  immunity  of  approximately  5 
dB-sm  demonstrates  that  classification  performance  can  be  significantly  improved. 
Figure  31  exhibits  the  smaller  noise  immunity  of  approximately  3dB-sm  associated 
Avith  channel  2  as  the  number  of  features  is  increased  from  3  to  10.  Figure  32 
demonstrates  an  almost  constant  increase  of  noise  immunity  of  3  dB-sm  to  5  dB- 
sm  as  the  number  of  features  is  increased  from  3  to  .5  to  10  to  20  in  the  channel 


with  an  unknown  multiplicative  component  (channel  3).  Figure  33  demonstrafe.s 
a  similar  increase  in  performance  for  the  unknown  additi^-e  component  channel. 

Figure.®  30  through  33  illustrate  that  classification  ijerfonnance  may  be  im¬ 
proved  by  increasing  the  number  of  feature.s  when  limited  to  a  2  MHz  frequency 


baiul,  irregardless  of  the  channel.  Channel  2  demonstrated  tlie  smallest  increase  in 
performance  of  the  three,  however,  the  increase  ij  rmoch  greater  than  that  cbseived 


by  Kamis  (2)  when  increasing  the  number  of  features  over  larger  frequency  bands 
of  8  MHz  and  22  MHz. 

Figures  34  and  35  show  tlie  rc  lationship  between  the  number  of  features  and 
classification  performance  over  4  MHz  and  8  MHz  frequency  bands  for  channel  3. 
Both  figures  indicate  an  additional  noise  immunity  of  approximately  3  to  5  dB 
as  the  number  of  features  is  increased  from  5  to  10.  A  smaller  increaise  in  per¬ 
formance  is  obtained  by  increasing  the  number  of  features  from  10  to  20  and  10 
to  15  as  shown  in  Figures  34  anc  35  respectively.  Therefore,  as  frequency  band 
is  increased,  channel  3  demonstrates  a  decrease  in  the  amount  .^f  improvement  in 
classificati.vn  performance  obtained  by  using  an  increeised  amount  of  features.  Note 
that  chanrscl  2  displays  this  same  property.  Figure  36  demonstrates  that  the  ad¬ 
ditive  component  channel  exhibits  the  same  properties  as»5lrannel  2  and  channel  3 

la? 

concerning  frequency  band,  number  of  samples,  and  classification  performance. 

Therefore,  Figures  30  through  36  indicaie  that  if  one  is  limited  to  a  small 
frequency  band,  it  may  be  worth  the  cost  of  the  additioned  processing  time  to 
obtain  ebliniates  at  10  or  20  frequencies  regardless  of  the  channel.  How-ever,  if  a 
larger  frequency  band  is  available,  the  additional  processing  cost  may  outweigh 
tl'.e  slight  increase  in  classification  performance  obtained  through  the  additional 
feature. . 
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Parameter  Under  Observation:  Number  of  Features 

CUiannel  Model:  Coherent  AWGN  Channel 

A5!:imuth  angle:  0 

Start  Frequency:  8.0  MHk 

Stop  Frequency:  10.0  MHz 

Number  of  Frequencies:  Varies 
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Figure  .30:  Coherent  chanuei  number  of  features  .study  (8.0-10.0  MHz). 


Parameter  Under  Observation:  Nu;nber  of  features 

Cliannel  Model:  Non- Coherent  AWGN  Channel 

Azimuth  angle:  0 

Start  Frequency:  8.0  MHz 

Stop  Frecpiency:  10.0  MHz 

Number  of  Frequencies:  V'aries 
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Figure  31:  Nor  Coherent  channel  number  of  features  studv  (.$.0-10.0  MHz). 
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Parameter  Under  01)servatioii:  Nuiiil»er  of  Features 

Channel  Model:  Multiplicative  Component  AWGN  Channel 

Azimuth  angle:  0 

Start  Freqnency:  8.0  MHz 

Stop  Freqiiency:  16.0  MHz 

Number  o*"  Frequencies:  Varies 


'.lumber  of  fea*ures  study 


5.3  Comparison  of  AWGN  Channels  Performance 


A  direct  conipr.risoii  of  llu*  four  '■hanaels  has  been  made  in  terms  of  clas¬ 
sification  perform?.:!ce.  Figure  -»7  demonstrates  tlie  dependence  of  classification 
performance  on  the  channel  given  a  2  Mliz  freqnenc.v  band.  A  significant  perfor¬ 
mance  degradrtion  of  at  leeist  lb  oBsin  is  expenenced  for  channel  2  with  respect  to 
channel  I.  The  classification  performance  o!  channels  -3  and  -*  is  further  degraded 
as  a  result  of  target  information  lost  in  the  normalization  processes  used  to  remove 
the  unknov.-n  t multiplicative  or  additive)  components. 

Figure  3<s  shows  the  resulting  performance  of  the  same  channels  over  a  wider 
-I  MKz  frequency  band.  The  performanre  of  channels  2.  3,  and  4  improves  sig¬ 
nificantly  with  respect  to  the  performance  of  channel  1.  .At  lower  noise  powers, 
channel  3  is  found  to  outperform  channel  4  by  a  surprisingly  wide  margin  of  ap- 
pro.ximately  2  dBsm.  Figure  39  compares  the  performance  of  the  4  ciiannels  over 
an  8  MHz  bandwidth.  The  additive  component  channel  now  outperforms  the  mul¬ 
tiplicative  channel.  Note  that  the  performance  of  channels  2,  3,  and  4  all  fall  within 
b  dBsm  of  each  other. 


Figures  37.  38,  and  39  together  indicate  that  channels  with  unknown  mul¬ 
tiplicative  or  additive  components  which  must  be  removed  will  always  degrade 
performance.  This  is  understandable,  since  information  concerning  target  size  is 
lost  in  the  nornicJization  process  used  to  remove  the  unknown  (multiplicative  or 
additive)  component. More  importantly,  these  figures  indicate  that  the  performance 
of  a  channel  with  an  unknown  additive  coniponert  (channel  4;  degrades  at  a  faster 
rate  with  a  narrowing  frequency  band  than  a  channel  with  an  unknown  multi¬ 
plicative  coinponcuf  (chaime!  3).  Figures  37  and  38  clearly  indicate  that  knowing 
the  diff“rence  in  ab.solntc  RC'b  between  frequencies  (channel  4)  will  not  provide 


T 


any  better  performance  over  knowing  oniy  the  relative  chaiige  in  RCS  (channel  3). 
Finally,  the  three  figures  indicate  that  without  phase  information  classification 
performance  decreases  dramat.cally  with  a  decrease  in  frequency  band. 

Figure:  40  and  41  show  t!ic  classification  performance  of  the  four  channels 
using  -S  and  1.5  features  respectively  in  the  S  MHz  frequency  band.  Again,  chan¬ 
nel  1,  which  utilizes  phase  infoniiation,  performs  significantly  better.  Note  that 
the  relative  performance  of  each  channel  with  respect  to  the  other  three  remains 
almost  constant  in  both  figures.  Therefore  ail  three  <'hannels  tend  to  suffer  the 
same  loss  of  performance  with  the  lowering  of  the  number  of  frequencies  in  this 
large  8  MRz  frequency  band. 

Figures  42  displa}'s  the  resulting  poor  jierformance  of  channels  3  and  4  when 
the  number  of  features  is  limited  to  three  and  the  frequency  hand  (2  MHz)  is 
limited.  The  9  dB-sm  increase  in  noise  power  immunity  of  channel  2  with  respect 
to  channel  3  at  the  20%  misclassification  level  is  one  of  the  largest  changes  in 
jierformance  found  between  the  two  channels. 

Figure  43  demonstrates  the  relatively  large  improvement  in  classification  per¬ 
formance  of  channel  3  with  respect  to  channel  2,  especially  at  noise  power  levels 
of  less  than  h  dB-sm.  This  increase  in  relative  classification  performance  demon¬ 
strates  that  the  performance  of  channel  3  degrades  faster  than  the  performance 
of  channel  2  when  the  number  of  features  is  decreased.  This  same  property  was 
noted  in  the  number  of  features  study.  Finally,  Figi:re  44  demonstrates  an  even 
smaller  decrease  in  classification  performa.nce  for  channel  3,  further  backing  the 
performance  loss  —  number  of  frequency  relationship  for  channel  3  at  small  band- 
widths. 

Figures  45  and  46  compare  the  channels  over  a  different  band  from  20  MHz  to 
24  MRz  and  20  MHz  to  28  MHz.  The  channels  demonstrate  the  .same  characteri.s- 


I 


?2 


It* 


i 


ins 


I 


h 


K 


m  W.^nMrMXAKmJL'^  jtT  A*  AT JfcT  AT MJlATLArJIAlirJ^rJW'J^r  MIT^  .-fc  JH  •"« *-»^TL*r« 


,TL*r«irPinui 


tics  over  these  bands  as  over  the  8  MHz  to  12  MHz  and  8  MHz  to  16  MHz  bands. 
Therefore  the  channel  characteristics  which  have  been  noted  do  not  appear  to  be 
dependent  on  the  exact  band  in  the  HF  spectrum  used,  but  only  depend  on  the 
size  of  the  band. 
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Paramet.pr  Under  Observation:  Channel 

Channel  Model:  Varies 

Azimuth  angle:  45 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  10.0  MHz 

Number  of  Frequencies:  10 


iO  tJ)  60 

NOSE  POWER  (DB-SM) 

Coherent  CUiannel 
Non- Coherent  Clhannel 
Multiplicative  Component  Channel 
Additive  Component  Channel 


Figure  .37:  Channel  study  over  2  MHz  band. 


Parameter  Ur.der  Ohservation:  Channel 

Channel  Model:  Varies 

Azimuth  angle:  45 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  12.0  MHz 

Number  of  Frequencies:  *0 


^iOISE  POWER  (OB-SM) 


Coherent  Channel 
Non- Coherent  Channel 
Multiplicative  CVjMiponent  Channel 
Additive  Component  Channel 


Figure  -38:  C4iannel  study  over  4  MHz  band. 
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Parameter  Under  Observation:  Cliannel 

Channel  Model:  Varies 

Azimuth  angle:  45 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  16.0  MHz 

Number  of  Frequencies:  10 


_ _ Coherent  Clhannel 

—  -  -  -  Non-Coherent  Channel 

—  -  —  Multiplicative  Component  Channel 

—  —  —  Additive  Component  C^hannel 


Figure  .39:  Channel  study  over  8  MHz  band. 
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Parameter  Under  Observation:  CHiannel 

Channel  Model:  Varies 

Azimuth  angle:  0 

Start  Frequency:  8.0  MIIz 

Stop  Frequency:  16.0  MHz 

Number  of  Freciuencies:  .*1 


NOISE  POWER  (OB-SM) 

_ Coherent  CHiannel 

-  -  .  -  Non-C-oherent  CUiannel 

—  -  —  Multiplicative  t^omponenf  Cdiannel 


Figure  '10:  Channel  study  for  n  IVcfiuency.  ^  MHz  i)and. 
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Parameter  Under  Observation:  Channel 

Channel  Model:  Varies 

Azimuth  angle:  0 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  16.0  MHz 

Number  of  Frequencies:  1.5 


Pi)  RO 

NOISE  POWER  (OB-SV) 


_ Coliereut  Channel 

-  -  -  -  Non-Coherent  CUiannel 

—  -  —  Multiplicative  C’omponent  Channel 


Figure  41:  Channel  study  for  i.5  frequency,  8  MHz  basjd. 
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Parameter  Under  Observation:  Channel 

Channel  Model:  Key 

Azimuth  angle:  0 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  10.0  MHz 

Number  of  Frequencies:  5 


Coherent  C-'hannel 


-----  Non- Coherent  Cihannel 
—  -  —  Multiplicative  C^omponent  Channel 


Figure  4.3;  C^hannel  study  lor  .t  frequency,  2  MHz  band. 
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Parameter  Under  Observation:  Channel 
Channel  Model:  Kej' 

Azimuth  angle:  0 
Start  Frequency:  8.0  MHz 
Stop  Frequency:  10.0  MHz 
Number  of  Frequencies:  10 


_ Coherent  Channel 

-  '  -  -  Non-Coherent  C^hannei 

—  -  —  Multiplicative  C'oinponeut  Channel 


Figure  44:  Channel  study  for  10  frequency,  2  MHz  band. 
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Parameter  Under  Obserratioii:  Cliaiinel 

Channel  Model:  Key 

Azimuth  angle:  45 

Start  Frequency:  20.0  MHz 

Stop  Frequency:  24.0  MHz 

Number  of  Frequencies:  iO 


_ Coherent  Channel 

-  -  -  -  Non-Coherent  CUiannel 

—  -  —  Multiplicative  Clomponenl  Channel 


Figure  45:  Channel  study  for  upiJer  4  MHz  baiul 
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CHAPTER  VI 

RELATIVE  PHASE  CHANNEL  MODEL 

6.1  Feature  Impiementation 

The  development  of  the  IT  feature  in  Sectioi:  3.5  was  intended  to  provide 
information  related  to  the  intrinsic  pkeise  of  the  target  at  times  when  the  accmate 
estimation  of  absolute  path  length  is  difficult  to  achieve.  It  has  previously  been 
shown  that  given  certciin  channel  conditions  one  may  create  ( —  5 )  target  features 
from  the  difference  of  the  intrinsic  phases  at  adjacent  frequencies.  These  (A^  -  1) 
features,  coupled  with  N  estimates  of  the  RGS  of  the  target  at  each  frequency, 
provide  a  total  {2N  -  1)  features  of  the  target  which  may  be  used  for  classification 
purposes.  The  (A'^  —  1)  IT  features  contain  additional  information  abotu  the  target 
provided  they  are  not  too  severely  affected  by  noise.  If  used  correctly  in  a  clas- 
.‘.oiraiion  algorithm  along  with  the  amplitude  features,  these  IT  features  should 
improve  classification  performauce  over  that  of  using  the  amplitude  features  only. 

A  v'ariation  of  the  nearest -neighbor  algorithm  was  developed  to  incorporate 
the  W  features  in  the  classification  process.  Previous  experience  has  shown  the 
nearest-neighbor  algorithms  can  provide  powerful  classification  performance  results 
(2),  and  therefore  have  been  used  as  the  starting  point.  In  addition,  the  use  of  a 
distance  me*ric  makes  intuitive  sense  in  a  problem  of  this  nn  ure.  Consider  two 
catalog  targets  [k  =  i,  2)  whose  backscatter  amplitudes  ( \/RCS)  measured  at  /i, 
/o,  and  /3  are  plotted  in  3  dimensional  space  as  in  Figure  47.  These  two  targets 
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have  two  corresponding  intritisic  M'’  values  each  given  by  Eq.  (6.1)  and  Eq.  (6.2) 
which  may  be  plotted  in  2-diniension<il  space  as  shown  in  Figure  48. 


W2}  -  rWih 

k  =  1,2 

(6.1) 

h 

j2 

11 

(6-2) 

where  %(/,  5  intrinsic  phase  of  target  K  at  frequency  f.. 

Ne.\t  cfsnsifier  e^ti mating  tl»e  r«>rr«j|yl«d  amplitude  and  corrupted  11  values  at 
an  unknown  target  u  at  the  corresponding  frequencies  /].  f2,  and  and  plotting 
the  estimates  in  Figure  47  and  Figure  48  as  shown.  Use  of  a  Euclidian  distance 
metric  provides  four  distinct  measures  as  given  in  Eq.  (6..3),  Eq.  (6.4),  Eq.  (('.5), 
and  Eq.  (6.6). 


<1  = 


<2  “  JLiMfi)  -  Au(fi))- 

1=1 


i=l 


<2  = 

j=l 


A  final  distance  between  the  unknown  target  and  the  two  catalog  targets  may 
then  be  constructed  as: 


(meters) 


(6.7) 


jA+W  _  jA  , 

®u,2  “  *u,2  +  PVOu.2 


(meters) 


where  the  p  term  provides  a  weighting  of  the  W  metric,  while  the  Q  term  provides 


a  unit  conversion  from  an  angle  measure  to  a  length  measure. 


The  Q  term,  labeled  the  cattilog  normalization  factor,  is  a  meastire  of  the 


separation  of  the  targets  in  the  catalog  set.  In  Figure  47,  the  average  separe-tion 


of  the  two  catalog  targets  in  terms  of  amplitude  features  may  be  calculated  by 


Eq.  (6.9). 


Dc=  ^MSi)-Ai(fi)f 


(meters) 


Likewise,  the  average  separation  of  the  catalog  targets  in  terms  of  \V  features  may 


be  calculated  by  Eq.  (6.10). 


.«=1 


(degrees) 


(6.10) 


The  ratio  of  these  two  terms,  as  given  by  Eq.  (6.11),  provides  the  catalog 


normalization  factor. 


^  L>C 

Uc 


(6.11) 


The  catalog  normalization  factor  may  be  viewed  in  the  following  manner.  Fig¬ 


ures  47  and  46  show  the  same  two  catalog  targets  in  terms  of  two  different  met¬ 


rics  sepsurated  by  two  different  averse  distances  as  calculated  in  Eq.  (6.9)  and 


Eq.  (6.10).  By  multiplying  the  IF  features  of  the  catalog  set  by  Q  (rescaling  the 


phase  feature  space),  one  can  find  the  new  average  separation  of  the  catalog  targets 


in  terms  of  the  IF  features  as: 


(6.12) 


#=  C?"E(H2(/,)-H'i(/i))= 


(6.13) 


Dll'  ^Q-  pfl' 


(6.14) 


Dll'  = 


(6.15) 


Therefore  multiplying  the  catalog  IT  features  by  Q  has  the  effect  of  convert¬ 
ing  the  IT  feature  space  to  an  amplitude  feature  space  by  using  average  target 
separation  as  the  conversion  factor.  Likewise,  the  rescaling  of  by  Q  as  shown 
in  Eq.  (6.7)  and  Eq.  (6.8)  has  the  same  resulting  effect. 

The  value  of  p,  the  IT  feature  weight,  reflects  the  emphasis  placed  on  the  IT 
features  in  the  classification  process.  For  p=0,  the  IT  features  are  ignored,  and  the 
classification  decision  is  solely  determined  by  the  amplitude  distance.  For  p  >1, 
the  normalized  IT  features  are  more  heavily  weighted  than  the  amplitude  features. 

In  order  to  determine  a  proper  value  of  p,  one  must  consider  the  radar  system 
from  which  the  estimates  of  amplitude  and  phase  are  made.  In  general,  the  de¬ 
gree  of  corruption  of  the  IT  features,  with  respect  to  the  amplitude  features,  will 
determine  the  emphasis  or  de-emphasis  of  the  normalized  IT-distance  measure. 

6.2  Channel  Simulation 


The  relative  phase  channel  model  was  simulated  as  follows.  The  IT  values 
of  the  catalog  set  were  first  created  using  the  complex  backscattered  coefficients 
previously  loaded  in  the  catalog  set.  The  normalization  factor  of  the  catalog  set 
was  then  created  by: 


^cnia/og  ~  ] 


(6.J6) 
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where  Aa{fi),  are  the  backscatter  amplitudes  of  catalog  targets  a  and  6, 

and  Wa(f,),  are  the  W  values  of  catalog  targets  a  and  b.  A  corrupt  test  set 

vector  was  formed  through  the  addition  of  coherent  zero  mean  Gaussian  noise  to 
a  high  accuracy  test  set  vector  as  described  earlier.  The  (N  -  1)  noisy  IT  features 
were  then  created  from  these  coherent  backscatter  amplitudes  by: 


=  fan' 


RM+l) 


'^•+1  yRrifi) 


i  =  l.iV  ~  1(6.17) 


where 


=  XTifi)  +  iV',  and 

RT(f,)  =  RTifi)  +  X", 

given  i?7'(/,),  Xj'{f,)  are  the  high  accuracy  complex  backscatter  coefficients  of 
te.st  target  T  at  frequency  /,•  (  or  /i+i),  and  N' ,  N"  are  i.i.d.  Gaussian  noise 
samples.  Note  that  the  distribution  if  TT'u(/J  is  going  to  b'"  denendent  on  Ax(fi) 
and  the  distribution  of  Au{fi).  This  may  be  visualized  in  Figure  49.  This  graphic 
description  of  the  addition  of  coherent  noise  to  a  large  and  small  target  shows  the 
strong  dependence  of  standard  deviation  of  on  the  uncorrupted  amplitude 

Axift)  of  the  data  point.  Therefore  targets  with  large  RCS  values  are  biased 
toward  having  W  features  that  are  less  corrupted  for  a  given  noise  pov/er. 

An  ujiderstanding of  the  relationships  between  the  distributions  of  W^ifi)  and 
Au{f,)  for  given  values  of  A^ifi}  and  levels  of  noise  power  is  important  in  order 
to  determine  sensible  values  of  p.  To  accomplish  this,  a  study  was  conducted  by 
first  selecting  targets  with  a  wide  range  of  backscatter  amplitudes.  The  standard 
deviation  of  the  IF  features  were  estimated  as  described  in  Degroot  [lOj  by: 

■  100 


(6.1$) 


where  H  j'(/j)  is  the  high  arcuracy  W"  feature  and  W^{fi)  is  the  rfh  independent 
corrupt  ir  feature.  The  min  statement  specifies  the  minimum  angle  between  the 
values  to  eliminate  branch  cut  effects.  In  a  similar  manner,  the  standard  deviation 
of  the  corrupted  amplitude  features  was  calculated  by  Eq.  (6.19). 


1 


""Mfi)  -  99 


100 


EiMfi)  -  AW>)f 


r=l 


(6.19) 


Estimates  of  j  versus  -4(/,  )  at  noise  power  levels  of  10  dB-sm  and  20  clB- 
sm  are  shown  in  Figure  -50.  Note  that  is  highly  dependent  on  both  .4(/J 

and  the  noise  power.  As  a  result,  one  may  desire  to  weight  the  distance  between  an 
unknown  IT  feature  at  frequency  /„  and  the  catalog  IT  feature  at  fn  dependent 
on  the  amplitude  feature  A|,  at  that  frequency.  The  standard  deviation  of  the 


amplitude  feature  (cr  i  )  was  found  to  be  almost  independent  of  A(/,)  and  only 


dependent  on  noise  power.  Estimated  values  were  >r^  ^  ^  0.3  sm  and  1.0  sm  at 


noise  power  values  of  10  dBsm  and  20  dBsm,  respectively.  The  previous  rtudy 
led  to  the  development  of  two  IT  feature  distance  measurements.  The  standard 
amplitude  distance  measurement  given  by; 


J.4  _ 


N 


Z{Ah(h)-An(L)y 

1=1 


(6.20) 


was  used  to  measure  the  distance  betw'een  catalog  target  k  and  unknown  target  it 
in  terms  of  amplitude  features.  One  IT  distance  measurement  was  given  by: 


1.  = 


’u,k 


N-1 


:=1 


(6.21 : 


where  n7.(/, )  and  Il'i,{/j)  represent  the  catalog  and  corrupted  IT  features  respec¬ 


tively.  .4  second  IT  distance  measurement  was  given  by: 

1 

fA-l  ..1 


A  -  I  \  /?(.4(..^,))(H '■(/,) -H'»(/,))' 

L:=i 


(0.22! 
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where  R{  .4(  /, ))  is  a  j'urictioix  of  the  corrupted  amplitude  of  the  test  vector.  R{  .4(  f, ) ) 
is  a  weight  which  approaches  one  for  large  values  of  -4(/, )  (<7  j, , ,  small)  and  ap- 
proacltes  0  for  small  values  of  .4(/, )  (err,,,  large).  Note  that  the  amplitudes  of 

-‘'v/i ; 

tlie  coriui>ted  vector  must  be  u.sed  and  not  the  amplitudes  of  the  catalog  pre.sently 
being  tested  against.  Use  of  the  uncorrupted  catalog  amplitudes  would  bias 
to  be  small  for  small  targets  and  large  for  large  targets,  therefore  biasing  the 
classification  decision. 

6.3  Experimental  Results 


(4assificatioii  runs  were  made  using  Eq.  (6.21)  for  to  determine  the 

effects  of  p  on  the  classification  algorithm.  Figures  .'>1.  -52,  and  53  show  the  results 
of  nine  tests  run  over  three  separate  frequency  bands.  Figure  .51  demonstrates 
an  average  improvement  of  b  dB-sm  of  noise  power  immunity  using  the  weight  of 
p  =  1.0  over  the  small  2  MHz  band.  However,  Figures  52  and  53  demonstrate  a 
loss  in  performance  as  the  j^.  distance  weight  is  increased.  These  results  indicate 
tha!  the  idj';.  distance  measure  is  less  than  optimal  since  it  weights  all  IF  features 
equally,  irregardless  of  the  amount  of  corruption  of  the  IF  feature. 

A  second  set  of  tests  were  performed  to  investigate  possible  functions  for 
/?(,4(/, ))  given  in  Eq.  (6.22).  In  this  set  of  tests  the  weighting  was  set  at  p  =  0.5. 
The  two  functions  chosen  fer  R[A{fi))  were: 

«(■•!(/;))  =  (0.23) 


and 


mMf,)) 


Aif,) 

.36.0 


(6.24) 
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since  a  value  of  .4(/,)  =  36  meters  corresponds  to  a  RCS  of  approximately  1300 
scjuar"  iiK'iers.  or  31  dB-sm.  corresponding  approximately  to  the  largest  RCS  value 
in  the  t  target  library  at  HF  frequencies. 

Figure  54  demonstrates  improvement  in  clcissification  performance  for  the  rel¬ 
ative  phase  cl  annei  over  the  non-coherent  channel,  however  the  increase  in  perh.M- 
niance  is  less  than  that  shown  in  Figure  51.  Figures  55  and  -56  show  little  or  no 
iiiiprovement  in  classification  perforiiiance  of  the  relative  phase  channel  over  the 
non  coherent  channel  when  utilizing  the  ->d„  l  measure.  Therefore  our  choices  of 
RiA(j,])  are  far  from  ideal. 

Overall  this  study  has  indicated  that  for  small  frecjuency  bands,  the  relative 
pilose  cliaiinel  described  may  improve  performance  over  the  non-coheient  channel 
if  the  relative  ph.ase  estimates  are  used  as  described.  However,  at  larger  frequency 
Ijdiids,  the  use  of  relative  phase  information  proposed  appears  only  to  deteriorate 
performance.  This  may  be  due  to  several  factors.  First,  the  development  of  the 
use  of  Hie  relative  {ihase  features  has  been  based  on  intuitive  judgement,  and  not 
proven  estimation  theorj.  Therefore  the  use  of  the  phase  features  while  making 
logical  sense,  may  be  far  from  the  optimum  use  of  relative  phase  as  a  classification 
feature. 

.4n  alternate  method  of  utilizing  relative  phase  information  involves  estima¬ 
tion  of  the  unknown  targets  phase.  Given  a  set  of  coherent  backscatter  amplitudes 
one  may  use  the  Cross  Correlation  .Algorithm  to  provide  a  maximum  likelihooi4 
estimate  of  the  unknown  target  phase,  tliereby  minimizing  the  condition.d  ri.'-k  for 
a  uniform  cost  function.  The  phase  estimate  may  then  be  used  to  aid  in  clas.sify- 
ing  the  target.  Research  has  /.hown  that  use  of  the  Cross  Correlation  .Alg<>ritiim 
will  coii^istciiMy  improve  classification  jierformance  over  u.'-e  of  amplitude  only 
estimates  jH|. 


i 


II  is  critical  to  note  these  results  are  highly  dependent  on  the  noise  model 
used  in  the  study.  Given  a  different  noise  model  where  the  noise  corrupting  the 
amplitude  and  phase  features  is  independent,  different  results  would  be  expected. 
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I  Paraim‘t<*r  Ihidcr  Observafion;  p 

Channel  Model;  Relative  Phase  Channel 
Azimuth  angle:  45 
Start  Frequency:  8.0  MHz 
I  Slop  Frequency:  10.0  MHz 

I  Number  of  Frequencies:  10 

I  H'  Distance  Measure:  id^j. 
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Pnraniefer  Uiulor  Ohs^rvaiion:  p 

Chamsei  Model:  Relative  Phase  Channel 

Azimuth  angle:  45 

Start  Frequency:  20.0  MHz 

Stop  Frequency:  24.0  MHz 

Number  of  Frequencies:  10 

11'  Distance  Measure: 

p  :  Key 
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P  =  0.05 


Figure  52:  11'  feature  high  4  MHz  band  weighting  study 
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Paranioler  Under  Observation:  p 

Channel  Model:  Relative  Phase  Channel 

Azimuth  angle:  30 

Start  Frequenc.v:  8.0  MHz 

Stop  Frequency:  10-0  MHz 

Number  of  Frequencies:  5 

ll"  Distance  Measure:  id„ 
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Figure  .^3:  lU  feature  low  -l^^illz  band  weighting  study. 
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Paranif'tpr  Uudor  Observation:  Tnclivirhial  Distanre  Weiglitings 

Channel  Model:  Relative  Phase  Channel 

Azinuith  angle:  45 

Start.  Frequency:  20.0  MHz 

Stop  Frequency:  24.0  MHz 

Number  of  Frequencies:  10 

IF  Distance  Measure:  2d,(  ^ 

p  :  0.5 
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Figure  55:  IF  feature  high  4  MHz  band  distance  weight  study. 
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CHAPTER  VII 


FADING  CHANNEL  MODEL  STUDY 

7.1  Implementation 

Multipath  propagation  is  a  common  occurrence  in  ionospherically  propagated 
signals.  In  the  radar  problem,  multipath  phenomena  generally  results  in  the  re¬ 
ceived  signal  being  the  result  of  the  two  or  more  signals  adding  in  and  out  of  phase 
at  the  receive  antenna.  One  signal,  which  travels  the  expected  path,  may  be  con¬ 
sidered  deterministic  in  nature.  The  other  signals,  which  travel  undesired  paths, 
may  best  be  described  as  a  stochastic  process  since  the  paths  they  take  through 
the  ionosphere  may  be  considered  less  stable.  As  a  result,  the  propagation  path 
may  best  be  described  as  a  fading  channel  [27]. 

A  fading  mode)  representing  the  estimation  of  features  from  multipath  prop¬ 
agated  signals  as  shown  in  Eq.  (7.1)  was  implemented  in  the  RSSE  program  [28]. 

j  ^  1,  Ar(7.1) 

where  aAji  fi)e^^Tifi)  is  the  deterministic  component,  Aj'(ft)Z{fi)e^^z(^) 

denotes  the  faded  component,  and  a  is  the  weight  of  the  deterministic  component. 
The  parameter  A'f(f,)t^^^^A  is  the  high  accuracy  complex  backscatter  component 
of  target  T  at  frequency  /j,  Z{i)  is  the  random  amplitude  of  faded  component 
taken  from  a  correlated  or  uncorrelated  rancloiu  process,  and  0^(i)  is  the  random 
phase  of  faded  component  taken  from  a  correlated  or  uncorrelated  random  process. 


Construction  of  the  simulated  feature  vector  is  shown  in  Figure  57.  The 

complex  backscatterecl  data  are  simply  scaled  by  the  chosen  variable  q,  where 

0  <.  a  <  1,  to  construct  the  deterministic  component.  Construction  of  the  faded 

components  are  more  coniplex.  The  magnitude  of  the  high  accuracy  backscatterecl 

•1  1 

data  is  first  scaled  by  (1  —  a:‘')2  to  ensure  the  vector  resulting  from  the  sum  of 
the  deterministic  and  faded  components  has  a  total  power  ecjual  to  the  test  set 
vector  used.  Next,  the  random  samples  Z{i)  are  constructed  by  first  generating 
samples  from  a  white  zero-mean  Caussian  random  process  with  a  variance  (i7“)  of 
one.  The.se  samples  are  then  passed  through  a  filter  with  an  impulse  response  h(i) 
where  i  =0  to  jV  ~  1.  The  output  process  has  an  autocorrelation  function  given 


Bzik)  =  h(i)h{k  -f  0,  -tV  +  1  <  <  iV  -  1  (7.2) 

i=Q 

as  described  in  [9j. 

The  total  power  in  the  process  {  Z(i)  }  is  then  given  by: 

N-l 

RzW=  (7.3) 

j=0 

Therefore,  by  selecting  filter  values  h(i)  such  that  the  total  power  of  the  process 
{  Z(t)  }  is  one,  the  resulting  power  of  the  faded  channel  vector  will  be  equivalent 
to  th.e  power  of  the  test  vector  used.  The  filter  values  used  were  constructed  from 
the  filter  function  by: 


/)(/)  = - -  (7.1) 

where  h(i)  is  the  filter  function,  and  ^  is  the  correlation  coefficient  selected  by  user 
the  so  that  /?2^(0)  equals  1  independent  of  the  value  .V  selected.  By  increasing 
'I'  from  zero  ({  Z{i)  }■  highly  correlated)  to  infinity  ({  Z{i)  }  nncorrelated)  the 
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upt’i  atoi  simulates  various  degrees  of  correlation  between  tlic  faded  components  of 
a  given  measurement  vector. 

The  faded  components  are  then  weighted  by  the  samples  from  the  random 
process  {  Z(i)  }.  Phase  is  assigned  to  each  random  faded  amplitude  by  taking 
sample.s  from  an  i.i.d.  process  whose  random  variables  are  uniformly  distributed 
over  the  inteival  -tt  to  ttt.  The  resulting  complex  deterministic  components  and 
complex  faded  components  are  tiien  summed,  producing  the  faded  channel  vec- 
t<jr.  Additive  white  Gaussian  noise  is  then  added  coherently  to  the  faded  channel 
vectoi  as  described  earlier.  The  simulated  vector  represents  the  outpu.t  of  a  faded 
(hanmi  with  additive  Gaussian  noise.  IBy  applying  the  techniques  described  in 
.sectiems  4.3.1  through  4.3.3  to  the  faded  channel  vector,  variations  of  the  faded 
channel  model  such  as  loss  of  pha.se  information  or  the  addition  of  an  unknown 
multiplicative  component  may  be  synthesized.  As  a  result,  three  faded  channel 
models  with  cliaracteristics  analogous  to  those  described  in  4.3.1  through  4.3.3 
and  titled  channel  IF,  channel  2F,  and  channel  3F  respectively  were  investigated. 

Three  studies  were  conducted  utilizing  th.e  faded  channel  model.  First  the  re- 
!ationshii>  between  the  deterministic  component  weight  (a)  and  the  classification 
performa.nce  for  the  three  channels  was  investigated.  Next  the  correlation  coeffi¬ 
cient  was  varied  while  maintaining  a  constant  deterministic  component  weight. 
Finally,  direct  comparison  between  the  three  faded  channel  models  was  conducted. 

7.2  Deterministic  Weighting  Study 

Figure  .‘jS  slunvs  the  the  degradation  of  classification  of  channel  Fl  as  is 
decreased.  Note  that  significant  degradatioi;  of  p."rfornianceis  only  found  at  o  -  O.'i. 
Figure  r)0  shows  a  significant  degradation  of  classification  peiformauce  at  o  0.73 
and  n  —  O.o  for  channel  F2.  Classification  usin^  channel  F3  demonstrates  an 
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approximate  10  percent  degradation  in  performance  as  o  is  decreased  from  1.0  to 
0.9,  and  from  0.9  to  0.75,  as  indicated  in  Figure  GO.  Thus,  these  figures  indicate  tiiat 
the  i)erformance  of  channel  F2  maj'  be  severely  degraded  by  a  low  deterministic 
weight. 

Figures  61,  62,  and  6-3  are  similar  to  the  previous  three  figures  with  the  ex¬ 
ception  of  a  smaller  frequency  band  of  2  MHz.  More  significant  degradation  of 
performance  for  channel  FI  is  noted  at  a=0.75  and  n— 0.5  as  noted  in  Figure  61, 
which  was  expected  with  the  smaller  frequency  band.  Figure  62  demonstrates  a 
decrease  in  performance  of  25  %  for  a  given  noise  power  as  a  is  varied  from  1.0 
to  0.9.  Likewise,  Figure  63  shows  the  same  large  decrease  in  the  ability  using 
channel  F3  as  the  deterministic  weight  is  decreased. 

Overall,  the  performance  tends  to  degrade  by  a  much  larger  amount  for  the 
same  changes  in  a  at  the  2  MHz  band.  Therefore  the  deterministic  weight  becomes 
an  increasingly  important  characteristic  of  the  faded  channel  when  the  frequency 
l)and  decreases. 

.4  third  set  of  Figures  (64,  65,  and  66)  were  constructed  using  one-half  the 
number  of  data  points  over  the  same  8.0  MHz  to  10.0  MHz  band.  Further  decreases 
in  classification  performance  with  decreased  deterministic  weighting  are  expressed 
for  all  three  channels  as  expected. 
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Paranietor  Piuler  Observation:  Dotenninistir  Weight. 

Chaimel  Model:  Coherent  Faded  Channel 

Azimuth  angle:  0 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  1“'  ''  MHz 

Number  of  Frequ-  ss:  10 

Determinislic  W(.  ■:  Key 

C'orrelation  Coeflit  nt:  0..5 


NOISE  POWER  (DB-SM) 


_ a  =  1.0 

- o  ==  0.9 

—  -  —  n  —  0.7.5 

—  —  —  a  =  0.5 


Figure  58:  Weighting  study  for  coherent  4  Mllz  band  faded  duinnel. 
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Parnnipfpv  Tinder  Observation:  Delermitiistic  Woigli* 

Channel  Model:  Non-Coherent  Faded  Channel 

Azimuth  angle:  0 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  12.0  MHz 

Number  of  Frequencies:  10 

Deterministic  Weight:  Key 

Correlation  Coefficient:  0.5 
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Parameter  Under  Observation:  Delerminislic  Weight 

Channel  Model:  Multiplicative  Component  Faded  Channel 

Azimuth  angle:  0 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  12.0  MHz 

Number  of  Frequencies:  10 

Deterministic  Weight:  Key 

Correlation  Coefficient:  0..5 


no  ClO 

NOSE  POWER  (DB-SM) 


_ o  =  1 .0 

- o  =  0.9 

—  -  —  (I  =  0.7.') 

—  —  —  o  =  O.-V) 


higure  60:  Weighting  study  for  multiplicative  1  MHz  band  faded  channel. 
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Paranipfrr  I  -Ifr  Ohservaiion:  Detorminisf ic  Wnglit 
Channel  Model:  Non-Golierent  Faded  Channel 
Azimuth  angle:  45 
Start  Frequency:  8.0  MHz 
Stop  Frequency:  10.0  MHz 
Number  of  Frequencies:  10 
Detenninistic  Weight:  Key 


Correlation  Coefficient:  0.5 


Parameter  Under  Observation;  Deterministic  Weight 

Cljannel  Model:  Multiplicative  C-omponent  Faded.  Channel 

Azimuth  angle:  45 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  10.0  MHz 

Number  of  Frequencies:  10 

Deterministic  Weight:  Key 

Correlation  Coefficient;  0.5 
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Paranirier  Under  Observation:  Deterniinislir  WeiKlit. 

Channel  Model;  Coherent  Faded  Channel 

Azimuth  angle:  45 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  10.0  MHz 

Number  of  Frequencies:  5 

Deterministic  Weight:  Key 

Correlation  Coefficient;  0.5 


N06E  POVKER  (OB-SM) 


a  --  1.0 
a  =  0.9 
Q  =  0.75 


Figure  64:  \\eighting  studj  for  coherent  5  feature  -  2  MHz  band  faded  channel. 
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7.3  Correlation  Study 


Figures  67,  68.  and  69  demonstrate  the  effect  of  the  correlation  of  the  faded 
components  on  the  performance  of  the  three  faded  channels  over  a  2  MHz  frequency 
Ijand.  In  this  Ijand,  channel  FI  does  not  demonstrate  a  decrease  in  classification 
performance  with  increasing  correlation.  Channel  F2  exhibits  a  4%  decrease  in 
performance  when  faded  components  become  highly  correlated  as  shown  in  Fig¬ 
ure  68.  Channel  F.3  again  shows  major  improvement  as  a  classification  feature 
when  the  correlation  increases  as  displayed  in  Figure  69. 

Figure  70  demonstrates  the  decrease  in  classification  performance  of  chan¬ 
nel  FI  resulting  from  an  increased  cori’elation  of  the  faded  components.  In  Fig¬ 
ure  71,  channel  F2  exhibits  a  significant  change  in  classification  performance  at 
^’=0.1,  however  the  overall  performance  is  quite  poor  in  all  cases  due  to  the  low 
weighting  (o  =  0.75)  of  the  deterministic  component.  Figure  72  shows  the  increase 
in  classification  performance  of  channel  -3  as  the  faded  components  become  more 
correlated.  The  significant  gain  of  .3  to  5  dB-sm  in  noise  power  immunity  indicates 
the  correlation  coefficient  may  be  an  important  parameter  in  a  faded  channel  w'hich 
exhibits  an  unknown  multiplicative  component. 

Figures  67  through  72  indicate  several  relations  between  classification  perfor¬ 
mance  and  correlation  of  faded  components.  Channel  Fl,  which  utilizes  intrinsic 
phase  information  tends  to  show  a  slight  decrease  in  performance  with  an  increase 
in  correlation.  This  may  possibly  due  to  the  corrupted  amplitude  of  the  target 
appearing  consistently  smaller  (or  larger)  at  each  frequency  due  to  the  correlation 
of  the  Z{i)  samples.  As  a  result,  the  target  may  appear  to  be  closer  to  a  smaller 
(or  larger)  target  in  terms  of  the  Nearest  Neighbor  .Algorithm,  and  therefore  be 
misclassified. 


The  tendency  of  channel  P2  to  increase  performance  as  the  faded  components 
become  uncorrelated  backs  up  the  previously  developed  theory.  However,  since 
channel  F2  removes  all  phase  information  and  passes  only  the  noisy  absolute  am¬ 
plitude  backscatter  components  to  the  classifier,  one  may  expect  to  see  a  greater 
decrease  in  performance  with  increasing  correlation  than  channel  1.  There  is  no 
strong  indication  of  this  behavior,  however. 

Cfiiannel  F.3  which  adds  the  unknown  multiplicative  component  which  must 
be  removed  as  described  in  section  4. .3. .3  demonstrates  a  significant  increase  in  per¬ 
formance  with  an  increase  in  the  correlation  of  the  faded  components.  As  reasoned 
earlier,  highly  correlated  faded  comjmnents  tend  to  make  the  target  look  larger  (or 
smaller)  across  the  band.  As  a  result,  the  average  RCS  of  the  unknown  target  will 
tend  to  be  larger  (or  smaller).  The  process  of  normalizing  the  corrupt  backscatter 
amplitudes  to  remove  the  unknown  multiplicative  component  also  tends  to  remove 
the  averayf  contribution  of  the  faded  component.  As  a  result,  the  normalized  am¬ 
plitude  feature  is  more  dependent  on  the  deterministic  component  of  the  faded 
channel,  resulting  in  better  classification  performance. 
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Paraiiif'lfr  Obs«*rvaficm;  Correlation 

Cliaiiiiei  Moilel:  Non-Coherent  Faded  Channel 

Aziniulh  angle:  0 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  12.0  MHz 

Nuiiiber  of  Frequencies:  10 

Deterministic  Weight:  0.7-5 

C^orrelation  Coefficient:  Key 
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Figure  71:  Correlation  .slvuly  for  non-coherent  4 
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Parameter  Under  Observation:  Correlation 

Channel  Model:  Multiplicative  Component  Faded  Channel 

Azimuth  angle:  0 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  12.0  MHz 

Number  of  Frequencies:  10 

Deterministic  Weight:  0.7.5 

Clorrelation  Coefficient:  Key 


_ =  0.1 
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Figure  72:  Correlation  study  for  multiplicative  1  MHz  band  faded  channel. 


7.4  Comparison  of  Faded  Channels  Performance 


A  niia!  bfudy  was  conducted  to  directly  compare  the  performance  of  the 
three  channels  for  different  correlation  coedicients.  detenninistit  weighis.  fre<juency 
bands,  and  feature  number  combinations.  Figures  7-3.  7-1.  and  77)  directly  compare 
the  three  feature  types  using  a  2  .\lHz  band  for  diflerent  values  of  a.  .Although  per¬ 
formance  is  relatively  poor  for  both  channels  Fl  and  F2.  channel  F2  significantly 
outijerforms  channel  Fl  for  o  -  0.7.v  and  o  -  0.5  (Figures  74  and  7-5  respectively). 
Channel  Fl  significantly  outperforms  the  other  two  for  all  values  of  o  as  shown. 

.A  comparison  of  the  three  feature  tvjres  for  various  values  of  'P  is  given  in 
Figures  7fi.  77.  and  7$.  Figures  76  and  77  demonstrate  that  for  a  small  (highly 
correlated  faded  components),  channel  F.3  may  actually  outperform  channel  F2. 
The  figures  also  demonstrate  the  strong  performance  of  channel  Fl  irregardless  of 
the  correlation  between  faded  components. 

Figures  79  and  SO  compare  the  three  faded  channels  for  two  diflerent  values  of 
t\.  Both  figures  demonstrate  the  superiority  of  channel  F'l  over  channel  F2  in  terms 
of  classification  performance.  In  fact,  the  performance  of  channel  F3  approaches 
that  of  channel  Fl  for  low  noise  power  in  Figure  80. 

Figures  <S1.  ?'2.  and  83  demonstrate  that  the  resulting  performance  from  using 
channel  F3  features  may  be  .significantly  greater  than  using  channel  F2  features 
irregardless  of  the  correlation  of  the  faded  components.  These  figures  also  reinforce 
the  strong  performance  of  channels  utilizing  pha.se  in.t^onnation.  It  is  important  to 
note,  however,  that  cla.-vsification  perfomunice  of  a  i  jiannol  with  an  unknown  mui 
tiplicative  component  (channel  F3)  may  approach  the  performan<e  of  a  ciianncl 
which  maintains  al)solute  amplitude  as  w»ll  as  target  phase  infoimation  tchaii- 


and  83  respectively,  the  performance  of  channel  F3  degrades,  but  still  provides 
relatively  good  classification  performance  at  low  noise  powers. 

Direct  comparison  of  the  three  faded  channels  in  term  of  classification  per¬ 
formance  resulted  in  several  interesting  results.  Most  importantly,  the  si.perior 
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performance  of  channel  FI  indicates  that  phase  information  concerning  the  target 
is  maintained  even  for  relatively  small  deterministic  weights  (q  =  0.5,0.75).  There¬ 
fore.  a  system  Avhich  can  accurately  estimate  intrinsic  target  phase  is  still  desired 
even  if  the  propagation  path  is  subject  to  a  degree  of  multipath  interference.  Sec¬ 
ond,  the  study  indicates  that  under  certain  frequency  band,  deterministic  weight, 
and  correlation  conditions,  the  removal  of  the  unknown  multiplicative  component 
from  the  output  of  channel  F3  provides  a  set  of  features  whose  classification  perfor¬ 
mance  rivals  that  of  channel  Fl,  and  is  better  than  the  performance  of  channel  F2. 
Frequency  bands  of  4  MHz,  deterministic  weights  of  0.75,  and  highly  correlated 
faded  components  tend  to  accentuate  the  normalized  features  compared  to  other 
conditions. 

It  is  important  to  note  that  given  a  propagation  path  modeled  by  these  con¬ 
ditions,  which,  are  far  from  ideal  in  the  RTI  sense,  classification  may  be  performed 
using  only  relative  RCS  values.  Furthermore,  having  the  ability  to  estimate  ab¬ 
solute  RCS  would  not  provide  an  increase  in  performance.  In  fact,  the  use  of 
the  corrupted  amplitude  estimates  without  normalization  would  seriously  degrade 
classification  performance  (Figure  81).  In  this  case,  even  the  ability  to  estimate 
target  phase  would  only  decrea.se  the  misclassification  percentage  by  four  or  five 
percent  at  low  noise  power  levels. 

Overall,  the  study  has  indicated  that  if  one  is  given  a  multipath  propagated 
.signal,  there  is  little  need  of  a  calibrated  reference  to  e.stimate  absolute  Rf'S. 
Rather,  a  stationary  uncalibrated  reference  may  be  utilized  to  estimate  relat've 
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Parameter  Under  Observation:  Channel 
Channel  Model:  Key 
Azimuth  angle:  45 
Start  Frequency:  S.O  MHz 
Stop  Frequency:  10.0  MHz 
Number  of  Frequencies:  10 
Delerministic  VVeight.:  0.75 
Correlation  Coefiicient:  0.5 
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Figure  74:  Parlially  faded  2  hUlz  band  channel  study. 


Parameter  Under  Observation:  Channel 
Channel  Model:  Key 
Azimuth  angle:  45 
Start  Frequency:  8.0  MHz 
Stop  Frequency:  10.0  MHz 
Number  of  Frequencies:  10 
Deterministic  Weight:  0.5 
Correlation  Coefficient:  0.5 
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Figure  75:  Highly  faded  2  MHz  band  channel  studv 
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Parameter  Under  01)servation:  Channel 

Channel  Model:  Key 

Azimuth  angle:  45 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  MHz 

Number  of  Frequencie.s:  10 

Deterministic  Weight:  0.9 

Correlation  Coefficient:  0.001 
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Pigui'e  (6:  Highly  correlated  2  MHz  baud  channel  study. 
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Parameter  Under  Observation:  Channel 
Channel  Model:  Key 
Azimuth  angle:  45 
Start  Frequency:  8.0  MHz 
Stop  Frequency:  10.0  MHz 
Number  of  Frequencies;  10 
Deterministic  Weight:  0.9 
Correlation  Coefficient;  oo 
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-  -  -  --  Non-Coherent  Faded  Channel 

—  -  —  Multiplicative  Component  Faded  Channel 


Figure  78:  Uncurrelated  2  MHz  baiul  channel  study. 
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Parameter  Under  Observa'iion:  Channel 
Channel  Model:  Key 
Azimuth  angle:  0 
Start  Frequency:  8.0  MHz 
Stop  Frequency:  12.0  MHz 
Number  of  Frequencies:  10 
Deterministic  Weight:  0.7-0 
Correlation  Coefficient:  0..5 
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Figure  79:  Partially  faded  -I  MHz  band  channel  study. 
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Parameter  Under  Observation:  Channel 
Channel  Model:  Key 
Azimuth  angle:  0 
Start  Frequency:  8.0  MHz 
Stop  Frequency:  12.0  MHz 
Number  of  Frequencies:  10 
Deierministic  Weight:  0.5 
C'orrelation  Coefficient:  0.5 
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Figure  80:  Highly  fadeil  •!  MHz  band  channel  studj 


Parameter  Under  Observation:  Channel 
Channel  Model:  Key 
Azimuth  angle:  0 
Start  Frequency:  8.0  MHz 
Stop  Frequency:  12.0  MHz 
Number  of  Frequencies:  10 
Deterministic  Weiglil:  0.75 
Correlation  C'oefficient:  0.1 
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Figure  Si:  Correlated  -1  MHz  band  rhannel  .studv. 


Parameter  Under  Oliservatiou:  CUiannel 
Channel  Model:  Key 
Azimuth  angle:  0 
Start  Frequency:  8.0  MHz 
Stop  Frequency:  12.0  MHz 
Number  of  Frequencies:  10 
Deterministic  Weight:  0.7.5 
Correlation  Coefficient:  2 
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CHAPTER  VIII 


CONCLUSIONS 
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A  study  of  existing  HF  radar  sj'sfems  with  respect  to  use  as  a  nieasureir.ent 
facility  has  been  condiicted.  Examination  of  a  HF  skywave  radar  system  has  led 
to  the  incorporation  of  two  classes  of  channel  models,  the  Gaussian  and  fading 
channel  models.  Further  variations  of  the  channels  within  each  class  simulated  the 
possible  restrictions  which  the  external  environment  and  propagation  path  might 
place  on  the  measurement  system  in  terms  of  the  estimation  of  various  target 
features.  As  a  result,  relationship  between  the  various  channels  and  classification 
performance  in  terms  of  various  radar  system  parameters  were  found. 

The  simulation  studies  have  demonstrated  that  classification  performance  is 
largely  dependent  on  the  size  of  the  frequency  band  for  all  channel  models  im¬ 
plemented.  Moreover,  the  results  have  indicated  that  when  tiie  frequency  band 
is  extremely  limited  (  i.  2  MHz),  classification  performance  may  be  improved  by 
increasing  the  number  of  frequencies  at  which  target  features  are  estimated.  Ther- 
forc  if  one  is  restricted  to  a  small  frequency*  band,  performance  may  be  increased 
at  the  cost  of  additional  processing  lime. 

The  fading  channel  models  simulated  the  elfect  of  muitipath  propagation  on 
the  RTl  system.  Results  indicated  that  without  phase  information,  the  e.--f  .mation 
of  absolute  RC'.S  feature  may  not  lead  tii  additional  classification  perforu.anre  em  r 
the  estimation  if  relative  RC.S  feature.^.  Therfiire.  tliere  would  be  no  advantage  in 
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having  a  calibrated  reference  which  would  allow  the  estimation  of  absolute  RCS- 
A  variety  of  channels  have  been  simulated.  In  order  to  further  evaluate  classi¬ 
fication  i)erformance.  data  from  existing  HF  radar  systems  is  desired.  Given  a  large 
set  of  Doppler  spectral  estimates  containing  familiar  aircraft,  several  techniques 
may  be  investigated  to  estimate  various  features  of  the  targets.  A  comparison  of 
these  e.stimates  would  provide  insight  into  the  statistical  properties  of  real  esti¬ 
mates.  In  addition,  a  comparison  of  these  estimates  with  those  simulated  from 
the  Elect roScience  Laboratory  data  base  would  provide  a  measure  of  how  accu¬ 
rately  the  backscatter  coefficients  measured  on  the  compact  range  represent  those 
measu~'*''  in  real  HF  :a<.tai  "y'item. 
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